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Abstract 


This  final  report  covers  the  results  of  investigations  extending  over  the  period  15  Jan 
2000  -  30  Nov  2002.  The  report  consists  of  two  parts,  in  Part  1  a  computer  model  of  a 
vircator  developed  at  the  Texas  Technical  University  is  developed  using  a  powerful  code 
called  MAGIC.  Part  1  ends  with  design  proposals  leading  to  an  increase  in  power  output 
of  52  %.  Part  2  is  primarily  concerned  with  the  influence  of  an  axial  magnetic  field  on  the 
performance  of  a  coaxial  vircator,  including  stability  investigations. 
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Summary 


This  is  the  Final  Report  covering  the  results  of  investigations  extending  over  an 
approximately  three  year  period,  15  Jan  99  —  30  Nov  2002.  The  Report  is  divided  into 
two  parts.  In  Part  1  computer  model  of  a  coaxial  vircator  is  discussed  in  some  detail, 
whereas  Part  2  is  solely  devoted  to  the  investigation  of  the  effect  that  an  axial  magnetic 
field  has  on  the  performance  of  such  a  vircator. 

In  Part  1  a  powerful  computer  code  called  MAGIC  is  used  to  model  a  coaxial  vircator 
under  development  at  the  Texas  Technical  University  (TTU)  in  Lubbock,  TX.  This 
required  close  collaboration  with  the  team  at  TTU  developing  the  device  and  involved 
exchange  visits.  The  first  seven  months  were  spent  on  establishing  the  contact  and 
developing  an  early  computer  model  of  the  device.  In  the  following  period  of  12  months 
the  model  was  refined  by  incorporating  ongoing  design  changes  introduced  at  TTU.  This 
direcdy  led  to  improved  accuracy  of  modelling,  the  frequency  of  oscillations  and  the 
output  power  now  being  reproduced  to  within  ±  5%.  The  remaining  period  of  16 
months  was  devoted  to  some  additional  design  alterations  and  a  detailed  investigation  of 
the  interaction  process  at  the  electron  cloud  level.  With  a  better  understanding  of  the 
operation  of  the  vircator  it  then  proved  possible  to  improve  the  power  output  of  the 
device  by  52  %  by  redesigning  the  sensitive  anode  area  of  the  tube. 

In  Part  2  the  effect  of  an  axial  magnetic  field  on  the  performance  of  an  axial  vircator  is 
discussed  in  considerable  detail,  the  results  of  the  investigations  confirming  the 
expectation  that  the  magnetic  field  hinders  the  performance  of  the  device.  In  the  process 
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careful  analysis  of  a  conducting  magnetic  diode  has  been  performed,  revealing  some 
hitherto  overlooked  solutions  resembling  those  obtained  in  the  absence  of  the  magnetic 
field. 

The  report  ends  with  suggestions  for  future  work  based  on  the  results  of  the  above 
investigations. 

It  should  be  added  that  the  above  investigations  generated  two  papers  (one  of  them  in  a 
special  issue  on  High  Power  Microwave  Generation)  and  five  contributions  at 
international  conferences. 
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Part  1. 


1.1  General  Introduction 

All  our  investigations  have  been  carried  out  in  close  cooperation  with  a  research  group  in 
the  Department  of  Electrical  Engineering  and  Physics  at  Texas  Technical  University  in 
Lubbock,  TX;  the  research  group  under  the  leadership  of  Prof.  Magne  Kristiansen  was 
then  in  the  process  of  developing  an  axial  vircator  [1,2]  and  our  task  was  to  provide  a 
suitable  computer  modelling  support,  primarily  by  using  MAGIC  [3].  Since  this  was  a 
cooperative  effort,  it  proved  important  to  establish  a  close  contact  with  our  colleagues  at 
TTU.  We  did  that  even  before  submitting  our  proposal.  Our  friends  in  Lubbock  helped 
us  greatly  by  sending  us  engineering  drawings  covering  the  design  of  their  vircator  tube 
and  we  in  turn  included  some  of  the  drawings  in  our  Proposal.  We  then  met  some 
members  of  the  group,  in  particular  Prof.  M.  Kristiansen,  at  the  ICOPS99  conference  in 
Monterey,  CA.  Since  the  award  of  the  grant,  we  have  further  developed  our  contracts  by 
one  of  us  (X.  Chen)  visiting  Lubbok  early  in  June  2000  after  attending  the  ICOPS  2000 
conference  in  New  Orleans,  LA.  The  visit  proved  to  be  most  valuable,  since  in  computer 
modelling  it  is  always  important  to  see  the  relevant  device  in  situ,  possibly  when  actually 
operating.  Dr.  Chen  was  warmly  welcomed  and  the  hosts  went  out  of  their  way  to  help 
and  make  the  visit  well  worthwhile.  In  turn  this  greatly  speeded  up  the  process  of  setting 
up  the  correct  computer  model  on  the  equipment  at  QMUL.  In  November  2001  we  had 
the  pleasure  of  a  visit  from  Prof.  M.  Kristiansen,  who  was  touring  Europe.  This  gave  us 
the  opportunity  of  introducing  the  staff  involved  in  computer  modelling  and  show  the 
substantial  computer  facilities  available  at  QMUL.  Finally  in  April  2002  Dr  Chen  and  Mr 
Toh,  a  post  graduate  student  chiefly  involved  in  the  vircator  modelling,  visited  TTU  on 
the  way  back  from  IVEC02  conference  in  Monterey,  CA. 
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Virtual  Cathode  Oscillators  (vircators)  have  been  studied  quite  extensively  over  the  last 
two  decades  [4-6].  In  principle  a  vircator  employs  an  intense,  mostly  relativistic  electron 
beam,  the  density  of  the  beam  greatly  exceeding  that  in  a  space-charge-limited  diode 
[7,8];  this  leads  to  the  formation  of  an  unstable  virtual  cathode,  which  then  generates 
relaxation  oscillations.  In  practice  the  device  is  capable  of  generating  pulsed  power  at 
levels  of  several  GW  and  in  the  GHz  range  of  frequencies.  The  main  merit  of  a  vircator 
is  its  simplicity  of  construction,  the  device  in  principle  comprising  a  cathode,  a 
semitransparent  anode,  usually  in  the  form  of  a  metal  mesh,  followed  by  an  output 
waveguide  and  a  window,  as  shown  schematically  in  Fig.  1,  no  magnetic  field  being 
required.  Since  a  vircator  uses  a  low  impedance  relativistic  electron  beam,  theoretically 
there  should  be  no  limit  to  the  power  level  at  which  it  can  operate,  nonetheless  its 
inherent  broadband  characteristic  and  a  wide  spread  of  electron  energies  renders  it  an 
oscillator  of  low  electronic  efficiency,  usually  below  5%.  Thus  the  principle  of  generating 
microwave  radiation  using  virtual  cathode  oscillations  is  conceptually  quite  simple, 
however  the  actual  process  is  highly  non-linear.  Consequently  all  the  earlier  attempts  at 
mathematically  analyzing  the  interaction  process  suffered  from  a  number  of  limitations 
[1,4,9-14].  It  would  appear  therefore  that  computer  modelling  presents  a  better  method 
for  the  study  and  understanding  of  vircators. 

As  a  first  step  in  our  modelling  it  seemed  prudent  to  agree  on  the  basic  geometry  of  the 
TTU  vircator.  Following  Fig.  4  of  ref.  1  and  some  additional  discussions  during  our  visit 
to  Lubbock,  the  basic  dimensions  of  the  device,  as  shown  in  Fig.  2,  have  been 
established.  Here  the  cathode  and  the  anode  are  respectively  in  the  form  of  two  coaxial 
and  partly  overlapping  cylinders.  This  constituted  an  important  development  of  a  simple 
axial  vircator  and  it  has  been  initially  suggested  some  time  ago  by  the  Russians  [17].  The 
advantages  of  such  a  system  are  twofold:  on  the  one  hand  the  area  where  the  virtual 
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cathode  is  likely  to  be  formed  is  better  defined  and  hence  easier  to  control,  on  the  other 
hand  the  generation  of  a  hollow  beam  removes  some  of  the  power  limitations  imposed 
by  the  existence  of  a  solid  beam.  The  dimensions  of  the  two  overlapping  cylinders  are  as 
follows:  the  outer  cylinder  of  130  mm  radius  acts  as  a  cathode  and  has  a  3  cm  wide  velvet 
band  on  the  inside  acting  as  an  emitter.  The  inner  cylinder  of  100  mm  radius  acts  as  an 
anode  and  is  made  from  a  10  pm  thick  aluminum  foil;  the  anode  then  continues  as  a 
cylindrical  output  waveguide  of  98  mm  internal  radius.  A  500  kV,  30  ns  Gaussian  pulse  is 
applied  to  the  diode  section  of  the  tube.  Microwave  oscillations  are  then  generated  near 
the  cathode/anode  section,  the  microwave  radiation  traveling  down  the  output 
waveguide  for  some  1600  mm.  The  waveguide  is  supported  at  one  end  by  a  vertical 
flange  and  at  the  other  end  it  has  a  plastic  window  which  is  permeable  to  microwave 
radiation.  In  addition  a  number  of  thin  axial  rods  are  placed  on  the  inside  of  the  anode 
drum  to  support  the  foil.  They  will  be  ignored  in  our  model  since  they  merely  affect  the 
degree  of  transparency  of  the  anode  and  form  an  integral  part  of  it. 
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1.2  Early  investigations  -  2D  model 

Early  investigations  cover  the  period  of  approx,  seven  months,  from  the  beginning  of  the 
grant  on  15  January  2000  till  31  July  2000.  In  our  original  Proposal  under  the  heading 
‘Statement  of  Objectives’  we  have  indicated  that  during  the  first  year  of  the  investigations 
we  intend  to  establish  close  co-operation  with  our  colleagues  at  TTU  and  using  MAGIC 
set  up  a  computer  model  of  the  vircator.  We  have  completed  the  initial  part  of  this  task 
in  the  first  seven  months  of  our  investigations  and  a  brief  summary  of  the  results  is  given 
below.  In  this  we  have  used  the  2D  version  of  MAGIC  for  the  sake  of  simplicity.  As 
pointed  out  in  Section  1.1  it  took  us  some  time  to  develop  the  computer  model  based  on 
Fig.  2,  MAGIC  being  quite  an  elaborate  code.  We  then  decided  to  perform  a  number  of 
simple  tests. 

1.2.1  Pulse  delay 

We  noticed  in  Fig.  6  of  ref.  1  that  there  was  a  time  delay  of  some  20  ns  between  the 
electron-beam  and  microwave  power  pulses,  the  latter  being  measured  with  a  probe  near 
the  output  window.  Using  MAGIC  we  have  decided  to  model  this  particular  aspect  of 
the  system  by  observing  the  propagation  of  a  wavepacket  down  a  circular  waveguide. 
Here  the  wavepacket  is  supposed  to  represent  a  microwave  pulse.  The  waveguide  was 
chosen  to  have  the  same  dimensions  as  the  cylindrical  output  waveguide  of  the  vircator. 
Also,  following  ref.  1  it  was  assumed  that  the  waveguide  operated  in  the  TE11  mode. 

A  Gaussian  2  GHz  microwave  pulse  was  then  applied  to  a  suitably  placed  dipole  at  one 
end  of  the  waveguide.  The  corresponding  transit  time  measured  at  the  other  end  was 
rmodei  =6.44  ns,  the  theoretical  value  for  a  TE11  mode  being  rtheoretlcai  =5.60  ns  (the 
actual  computations  appear  on  p.  6  of  ref.  18);  as  mentioned  above  the  experimental 
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value  of  the  delay  was  rexp  =  20  ns.  The  difference  between  the  two  sets  of  results  could 

be  explained  by  the  fact  that  in  practice  it  takes  approx.  5-10  ns  for  the  microwave  pulse 
to  be  formed  and  the  oscillations  to  begin. 

This  simple  measurement  acted  as  a  primitive  test  of  our  model  and  also  confirmed  our 
correct  use  of  MAGIC. 

1.2.2  Frequency  spectrum 

Having  gained  some  confidence  in  the  more  rudimentary  performance  of  our  model  we 
then  decided  to  investigate  the  actual  microwave  output  and  its  frequency  spectrum.  The 
radial  field  component  Ep  as  a  function  of  time  and  the  corresponding  FFT  are 

respectively  shown  in  Figs.  7  and  8  of  ref.  1,  the  measurements  being  taken  with  the  help 
of  a  probe  8  mm  long  situated  1.5  m  down  the  waveguide.  In  our  simulation  experiment 
we  have  managed  to  repeat  the  above  result  with  limited  degree  of  accuracy.  Since  the 
vertical  scale  in  Fig.  7  of  ref.  1  is  not  marked  we  could  not  compare  the  amplitude  of 
oscillations  with  our  Fig.  3a,  but  the  waveform  of  the  microwave  output  is  fairly  similar 
in  the  two  cases.  The  FFT  of  the  output  waveform  is  shown  in  Fig.  3b.  We  find  that  the 
largest  frequency  peak  occurs  at  2.7  GHz  as  against  2  GHz  in  Fig.  8  of  ref.  1.  The 
difference  of  35  %  is  almost  certainly  due  to  some  structural  simplifications  of  the  early 
model  and  possibly  due  to  the  use  of  a  2D  version  of  MAGIC.  As  we  shall  see  later  this 
difference  was  finally  reduced  to  2.5  %  when  we  have  used  a  more  detailed  model  and  a 
3D  version  of  MAGIC.  As  an  additional  precaution  we  moved  our  Ep  observation  point 

from  the  circumference  to  the  centre  of  the  waveguide,  but  the  results  were  similar  to 
those  shown  in  Figs.  3a  and  b. 
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1.2.3  Output  waveguide 

As  one  more  preliminary  test  of  our  model,  we  investigated  the  electromagnetic  field 
distribution  inside  the  output  waveguide.  We  noted  that  the  mode  pattern  remained  fairly 
coherent  only  at  the  beginning  of  the  pulse  and  it  disintegrated  quite  rapidly  afterwards. 
This  was  to  be  expected  for  two  reasons:  vircators  are  well  known  for  their  broad-band 
output,  Fig.  3b  being  a  good  example.  Secondly,  as  has  been  pointed  out  before  the 
output  waveguide  is  heavily  overmoded,  a  situation  which  encourages  the  transfer  of 
energy  to  higher  order  modes. 

We  find  from  Figs.  4a  and  5a  that  in  the  early  stages  of  oscillations  the  two  dominant 
modes  appear  to  be  respectively  TM01  and  TE01.  The  early  appearance  of  the  TM01 
mode  should  not  surprise  us  in  view  of  the  existence  of  a  strong  Ep  field  component 

near  the  anode  surface.  The  appearance  of  the  TE01  mode  could  be  caused  by  some 
asymmetry  in  the  virtual  cathode  region  giving  rise  to  azimuthal  components  of  the 
electric  field.  However  it  is  not  clear  how  it  could  be  generated  in  our  model,  since  at  this 
stage  we  have  been  using  a  2D  version  of  MAGIC.  However,  the  above  results  differed 
from  the  experimental  probe  measurements  at  TTU,  where  they  were  inclined  to  deduce 
the  presence  of  a  TE11  mode  at  the  end  of  the  waveguide.  The  test  also  confirmed  the 
versatility  of  MAGIC  as  a  computer  code  for  modelling  microwave  devices. 
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1.3  Intermediate  investigations 


1.3.1  Introduction 

Intermediate  investigations  cover  the  second  year  of  the  grant,  starting  on  1  August  2000 
and  ending  on  31  July  2002.  In  order  to  improve  the  convergence  of  our  results  with 
those  obtained  at  TTU  we  have  carefully  incorporated  in  our  model  all  the  structural 
changes  dictated  by  experience  and  introduced  by  our  colleagues  at  TTU  during  the  first 
seven  months  of  our  cooperation.  In  addition  we  gained  enough  confidence  in  using 
MAGIC  to  change  from  2D  to  3D.  As  we  shall  see  both  changes  had  quite  a  dramatic 
effect  on  the  accuracy  and  reliability  of  the  results. 

1.3.2  Improved  TTU  vircator  model  -  3D  model 
1.3.2a  Design  alternations 

A  number  of  alternations  in  the  design  of  the  vircator  have  been  introduced  at  TTU 
during  the  first  seven  months  of  our  investigations,  the  final  design  of  the  important  gun 
section  being,  shown  in  Fig.  6.  In  particular  a  rectangular  septum  has  been  added  in  order 
to  improve  the  performance  of  the  device.  The  corresponding  model  of  the  final  version 
of  the  device  is  shown  in  Fig.  7a,  the  septum  being  shown  separately  in  Fig.  7b.  The 
overall  system  now  consists  of  a  Marx  generator  followed  by  a  pulse  forming  line 
feeding  two  coaxial  electrodes  through  a  spark  gap.  The  vircator  itself  comprises  a 
coaxial  diode  in  the  form  of  a  metal  cylinder  supporting  a  velvet  emitter  of  132  mm 
inner  radius  and  an  internal  cylindrical  anode  made  from  a  stainless  steel  mesh  of  60 
%  transparency  and  of  98  mm  radius.  At  one  end  the  mesh  anode  is  strengthened  by  a 
metal  wire  ring  (wire  radius  2.5  mm)  and  at  the  other  end  it  is  attached  to  a  flange  and 
a  1 .5  m  long  cylindrical  waveguide  which  has  an  output  window  at  the  far  end.  Two 
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orthogonal  probes  are  placed  at  the  end  of  the  waveguide  for  both  mode  detection  and 
power  measurement. 

In  operation  a  single  pulse  of  500  kV  and  70  ns  duration  is  applied  across  the 
electrodes.  Typical  experimental  input  voltage  and  current  pulses  are  shown  in  Fig. 
22,  where  the  measured  peak  values  are  respectively  731  kV  and  36.7  kA.  An 
indication  of  the  input  impedance  is  obtained  by  dividing  peak  voltage  by  peak 
current  which  gives  19.9  Q.  Fig.  23  shows  typical  experimental  input  and  microwave 
output  power  pulses,  their  respective  peak  values  being  24.7  GW  and  558  MW.  The 
quality  of  our  computer  modelling  using  MAGIC  will  be  judged  by  the  degree  of 
convergence  with  the  above  parameters. 

1.3.2b  Computational  details 

Let  us  now  consider  some  computational  details  of  the  simulation  process.  The  MAGIC 
code  is  very  flexible  and  therefore  well  suited  to  achieve  the  required  accuracy.  Since  the 
vircator  exhibits  a  large  degree  of  axial  symmetry  it  was  natural  to  choose  polar  (r,(p,z)  co¬ 
ordinates  for  our  modelling.  Since  we  are  now  using  a  3D  version  of  MAGIC  the  total 
number  of  cells  has  gone  up  by  a  factor  of  80,  which  is  approx,  the  number  of  cells  in  the 
azimuthal  direction.  Although  the  size  of  8<p  and  8z  cells  is  fixed,  the  8r  cells  can  be  graded 
by  choosing  smaller  values  in  regions  of  special  interest.  The  choice  of  cell  size  governs  the 
maximum  permissible  time  step  which  for  reasons  of  stability  must  satisfy  the  COURANT 
criterion 


XT 
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Here  c  is  the  velocity  of  light  and  5t  the  time  step.  Since  the  total  number  of  cells  is  quite 
large,  the  overall  time  of  a  simulation  run  frequently  last  as  much  as  24  hours,  even  with 
powerful  modem  PC  computers. 

In  practice  the  emitter  is  in  the  form  of  a  bobin  of  velvet  glued  to  the  inside  of  the 
cathode  cylinder,  its  action  being  that  of  a  field  emitter.  The  EXPLOSIVE  EMITTER 
command  in  MAGIC  is  well  suited  to  simulate  the  performance  of  such  a  cathode.  In 
our  model  it  is  assumed  that  the  thickness  of  the  emitter  and  its  emission  efficiency  are 
both  uniform  around  the  circumference.  The  anode  of  the  vircator  is  in  the  form  of  a 
light  metal  mesh  supported  by  eighty  thin  metal  rods  (not  shown)  and  a  thick  metal  ring. 
The  simulation  of  the  anode  structure  is  based  on  another  MAGIC  command  called 
FOIL,  the  command  ensuring  60  %  transparency  to  the  electron  stream  in  close 
agreement  with  experimental  results. 

In  order  to  improve  the  performance  of  the  vircator  a  rectangular  septum  has  been 
introduced,  as  shown  in  Figs.  6,  7a  &  b.  In  principle  there  are  two  ways  of  modelling  the 
septum:  we  can  use  either  Cartesian,  Fig.  8a  or  polar  coordinates.  Fig.  8b.  When 
Cartesian  coordinates  are  used,  the  septum  has  smooth  edges  but  the  cylindrical  walls  of 
the  waveguide  are  quite  rough,  the  accuracy  of  modelling  now  strongly  depending  on  the 
size  of  8x  and  8y  cells.  In  polar  coordinates,  the  cylindrical  walls  remains  smooth,  but  the 
shape  of  the  septum  is  now  critically  dependent  on  the  size  of  8(p  and  8r  cells.  Since  we 
have  decided  to  use  polar  coordinates  anyhow,  it  was  only  natural  to  follow  Fig.  8b  in 
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modelling  the  septum.  Also  the  angular  cell  size  8(p  was  chosen  to  be  0.0785  radians, 
which  corresponds  to  80  cells  around  the  circumference  of  the  cylinder.  The  use  of  80 
cells  for  the  angular  resolution,  together  with  a  large  number  of  8r  and  8z  cells  required 
computer  runs  often  in  excess  of  12  hours. 

We  now  have  to  consider  the  problem  of  power  supplies  and  the  modelling  of  the  feed 
system.  At  TTU  a  Marx  generator  supplies  a  500kV,  70  ns  pulse.  In  our  simulations  we 
use  a  voltage  pulse  of  the  same  amplitude  but  with  an  exponential  rise  time  given  by 

V(t)  =  V0  (l  -exp(-t/T0)) ,  V,  =  500kv ,  T0  =  4.0  ns  (1.2) 

followed  by  an  abrupt  fall.  Figs.  9a-d  show  different  ways  of  connecting  power  supplies 
to  the  system.  The  first  three  have  been  rejected  in  favour  of  the  fourth,  since  they  are 
too  near  to  the  excitation  (virtual  cathode)  region  and  additionally  in  Fig.  9c,  the  feed  is 
orthogonal  to  the  desired  direction.  Accordingly  Fig.  9d  with  its  coaxial  feed  geometry 
and  clear  separation  from  the  general  excitation  region  is  the  preferred  method.  Here  the 
length  of  the  separate  feed  channel  needs  to  be  taken  into  account  when  comparing  the 
delay  between  the  excitation  and  the  output  pulse  appearing  at  the  end  of  the  waveguide. 

Lasdy  let  us  consider  TTU  measurement  techniques.  Two  perpendicular  probes  situated 
1.5m  down  the  wave  guide  were  used  to  monitor  the  output  with  the  help  of  a  Hewlett 
Packard  8GHz  digital  oscilloscope  HP54720D  and  an  HP8719C  network  analyser.  The 
microwave  power  generated  by  the  vircator  was  then  obtained  by  taking  into  account  the 
attenuation  of  the  probes,  the  cables,  the  connectors,  the  fixed  attenuators  and  the  crystal 
detectors.  The  corresponding  frequency  spectrum  was  obtained  by  applying  a  fast 
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Fourier  transform  (FFT)  to  the  waveform  directly  recorded  by  a  digital  oscilloscope.  As 
an  example  typical  experimental  results  obtained  at  TTU  are  shown  in  Figs.  10a  and  b. 

At  this  point  we  should  compare  Figs.  3b  and  10b,  the  difference  being  mostly  due  to  the 
introduction  of  the  rectangular  septum.  The  effect  is  quite  dramatic,  the  frequency 
spectrum  in  Fig.  10b  being  much  cleaner  than  in  Fig.  3b,  the  septum  presumably  acting 
as  a  quasi-cavity.  At  the  same  time  the  main  resonant  frequency  was  shifted  downwards 
from  2.7  to  1.834  GHz. 

1.3.3  Testing  of  the  new  model 

Having  obtained  information  on  the  latest  changes  in  the  design  of  the  vircator,  as 
explained  in  Section  1.3.2a,  we  have  decided  to  put  our  new  model  to  the  test  by 
repeating  the  latest  results  conveyed  to  us  by  the  TTU  team  and  respectively  shown  in 
Figs.  10a  &  b. 

1.3.3a  Power  and  efficiency 

Let  us  consider  a  particular  simulation  run,  using  the  feed  system  of  Fig.  9d  where  a 
voltage  of  500kV  would  be  applied  between  cathode  and  anode  resulting  in  an  input 
power  characteristic  shown  in  Fig.  11;  here  the  value  of  T0  in  (1.2)  was  chosen  to  be  4.0 
ns(rise  time  characteristic),  the  input  reaching  a  steady  state  input  power  of  7.80  GW.  In 
Fig.  12  we  have  the  microwave  output  power  measured  over  the  cross-section  of  the 
waveguide  near  the  window.  There  is  a  delay  of  some  20  ns  in  the  rise  of  output  power 
due  to  the  time  it  takes  for  the  pulse  to  travel  from  the  generator  through  the  coaxial 
feed  system,  the  emitter  anode  space  and  along  the  waveguide,  before  reaching  the 
observation  probes.  The  shape  of  the  output  power  pulse  roughly  follows  that  of  the 
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input  power.  Taking  the  value  of  160  MW  for  the  output  power  the  device  efficiency  is 
of  the  order  of  2.0%  —  a  typical  figure  for  a  vircator  and  very  close  to  the  experimental 
value  quoted  by  TTU  (see  also  Fig.  23).  The  resulting  radial  component  of  the  electric 
field  Ep  at  a  point  near  the  centre  of  the  exit  plane  of  the  waveguide  can  be  seen  in  Fig. 
13a,  its  frequency  spectrum  (FFT)  being  shown  in  Fig.  13b. 

The  above  results  allow  us  to  compare  the  experimental  and  simulation  performance  of 
the  device.  In  Figs.  10a  and  13a,  we  have  the  electric  field  component  Epat  the  output 
plane.  Allowing  for  the  fact  that  a  Gaussian-shaped  input  voltage  pulse  was  applied  in  the 
experimental  case  and  a  steady  voltage  reached  by  an  exponential  rise  in  the  simulation 
runs,  there  is  reasonable  agreement  between  the  two  results.  Figs  10b  and  13b 
respectively  show  frequency  spectra  of  the  two  cases.  It  can  be  seen  that  not  surprisingly, 
a  ‘cleaner’  response  is  obtained  in  computer  simulation  than  in  the  actual  experimental 
results;  also  the  main  resonant  frequencies  agree  remarkably  well  to  within  1%.  All  this 
indicates  that  good  accuracy  can  be  obtained  by  using  MAGIC  for  modelling  the 
operation  of  a  vircator. 


1.3.3b  Size  of  cathode  to  anode  gap 

In  this  brief  test  we  have  investigated  the  influence  that  the  size  of  the  cathode  to  anode 
gap  has  on  the  frequency  of  oscillations.  In  order  to  allow  for  the  thickness  of  the  velvet 
emitter  its  height  was  set  to  be  equal  to  two  radial  cells,  as  show  in  Figs.  9a  —  d.  However, 
this  effectively  reduces  the  cathode  to  anode  gap  by  ~  1.0  cm  and  causes  an  upward  shift 
in  the  oscillation  frequency  [19].  We  have  therefore  decided  to  investigate  the  effect,  the 
results  being  summarised  in  Table  1.1. 
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Table  1. 1  Effect  of  the  cathode  to  anode  gap  on  the  resonant frequency. 


Height  of  emitter 

1.0cm  (two  cells) 

0.5  cm  (one  cell) 

0.0  cm  (0  cell) 

Cathode/ anode  gap 

3.2  cm 

3.7  cm 

4.2  cm 

Main  frequency 

1.988  GHz 

1.963  GHz 

1.816  GHz 

One  can  see  that  the  operation  frequency  of  the  vircator  is  critically  dependent  on  the 
size  of  the  cathode  to  anode  gap.  When  the  gap  is  close  to  that  in  the  experimental  set 
up  (zero  height  of  the  emitter),  we  find  from  Figs.  10b  and  13b  that  the  simulated 
frequency  agrees  to  within  1%  with  the  measured  one. 


1.3.3c  Presence  of  the  anode  ring 

A  thick  metal  ring  is  mounted  at  the  far  end  of  the  support  rods  (not  shown)  to  add 
strength  to  the  fragile  anode  mesh,  as  shown  in  Figs.  6  and  7a.  When  the  ring  is  removed, 
there  is  a  slight  upward  shift  in  the  resonant  frequency,  as  shown  in  the  last  line  of  Table 
1.2. 


Table  1.2  Frequency  shift  when  the  anode  support  ring  has  been  removed. 


Height  of  emitter 

1.0  cm  (two  cells) 

0.5  cm  (one  cell) 

0.0  cm  (0  cell) 

Cathode/ anode  gap 

3.2  cm 

3.7  cm 

4.2  cm 

Main  frequency 

2.138  GHz 

2.045  GHz 

1.853  GHz 
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Again  we  find  that  the  removal  of  the  anode  support  ring  has  a  direct  bearing  on  the 
resonance  frequency  of  the  vircator.  The  results  of  Tables  1.1  and  1.2  are  plotted  in  Fig. 
14.  One  can  see  that  the  size  of  the  cathode/anode  gap  would  affect  the  transit  time  of 
the  electrons  and  thus  the  oscillation  frequency  of  the  vircator.  Similarly  the  presence  or 
absence  of  the  support  ring  would  strongly  affect  the  electric  field  in  the  sensitive 
interaction  region  of  the  vircator  and  thus  influence  the  electron  dynamics  of  the  system. 
Thus,  as  was  to  be  expected,  the  above  investigations  have  confirmed  the  fact  that  the 
operation  frequency  of  the  vircator  is  critically  influenced  by  the  fields  prevailing  in  the 
interaction  region  and  the  related  electron  dynamics  (see  also  Section  1.3.3d  and  1.4.1). 
This  point  further  emphasizes  the  difficulties  associated  with  attempts  to  investigate  the 
performance  of  a  vircator  analytically  (see  also  Section  1.4.2d). 


1.3.3d  Position  and  size  of  the  septum 

As  shown  in  Figs.  6  and  7a,  b  a  rectangular  septum  has  been  added  to  the  TTU  vircator. 
In  general  a  septum  can  be  regarded  as  a  circuit  element  designed  to  improve  the 
performance  of  a  vircator  by  forming  a  quasi  cavity.  In  Fig.  6  the  distance  between  the 
vertical  part  of  the  cathode  structure  and  the  septum  is  11.4  cm;  this  would  constitute  a 
half-wavelength  at  1.59  GHz  (TE11  mode)  or  at  1.76  GHz  (TM01  mode),  the  latter 
being  nearer  to  both  measured  and  simulated  values.  There  are  however  other  factors  to 
be  taken  into  consideration,  principally  the  presence  80  very  thin  anode  supporting  rods 
(not  included  in  the  model)  and  the  electron  cloud  which  would  alter  the  resonant 
frequency  by  introducing  a  reactive  term. 

A  number  of  computer  runs  have  been  made  with  the  septum  at  different  positions 
along  the  anode,  the  results  being  summarised  in  Table  1.3  and  Fig.  15. 
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Table  1.3  Effect  of  septum  position  on  the  resonant frequency. 


Location  of 

the  septum 

4.6  cm 

5.6cm 

6.6  cm 

*7.6  cm 

8.6  cm 

‘cavity’ 

length 

14.4cm 

13.4  cm 

12.4  cm 

11.4  cm 

10.4  cm 

Operation 

Frequency 

1.755  GHz 

1.772  GHz 

1.806  GHz 

1.816  GHz 

1.858  GHz 

*  The  original  location  in  TTU  experiments 


The  shift  in  frequency  is  as  we  would  expect  —  a  decrease  in  frequency  is  associated  with 
an  increase  in  the  length  of  the  ‘cavity’  formed  between  the  cathode  base  and  the  septum. 
This  simple  test  confirms  that  in  the  presence  of  a  septum  the  resonant  frequency  of  the 
vircator  critically  depends  on  its  position  along  the  axis  of  the  tube. 

Another  septum  featuring  a  larger  area  and  a  smooth  edge,  as  shown  in  Fig.  1 6,  has  also 
been  examined.  A  brief  comparison  of  the  results  is  given  in  Table  1.4. 


Table  1.4  Effect  of  the  septum  siye  on  the  operation  of  the  vircator. 


Septum  size 

Original 

Large 

Spectrum 

1.816  GHz 

1.882  GHz 

Efficiency 

2.05  % 

1.54% 

The  corresponding  frequency  spectrum  with  the  new  septum  in  place  is  shown  in  Fig.  17; 
the  output  is  rather  disappointing,  showing  a  pronounced  peak  at  1.667  GHz  in  addition 
to  the  main  resonant  frequency  at  1.882  GHz.  Also,  the  efficiency  has  fallen  from  the 
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usual  value  of  2.05  %  to  1.54  %,  probably  due  to  excessive  blocking  of  the  RF  power 
transmitted  down  the  waveguide  (impedance  mismatch). 
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1.4  Final  Investigations  -  3D  model 

This  period  covets  the  last  16  months  of  our  investigations,  that  is  from  1st  August  2001 
till  30  Nov  2002.  We  have  now  reached  a  stage  when  our  confidence  in  the  reliability  of 
our  model  was  sufficient  to  attempt  an  overall  simulation  of  the  performance  of  the  TTU 
vircator 


1.4.1  Size  and  position  of  circular  septum 
1.4.1a  General  comments 

During  our  visit  to  TTU  at  the  end  of  April  2002  we  have  been  asked  to  investigate  the 
effect  of  a  septum  in  the  form  of  a  circular  disc  with  a  rectangular  cut-out,  on  the  general 
performance  of  the  vircator.  This  request  was  very  much  in  line  with  our  overall  desire  to 
bring  investigations  based  on  computer  modelling  and  MAGIC  as  close  as  possible  to 
the  design  requirements  of  the  device  and  thus  to  achieve  the  main  objective  of  being 
able  to  treat  computer  modelling  as  a  genuine  design  tool.  This  can  only  be  achieved 
when  the  code,  in  our  case  MAGIC,  reaches  a  sufficient  degree  of  sophistication  and 
reliability  to  be  able  to  replace  machine  modelling. 


1.4.1b  Dimensions  and  the  simulation  setup 

An  axial  cross-section  of  the  interaction  region  of  the  TTU  vircator  with  the  new  septum 
is  three  different  positions  are  shown  in  Fig.  18a  where  the  dimensions  (in  mm)  and 
relative  positions  of  the  cathode  cylinder,  the  electron-emitting  strip  of  velvet  and  the 
anode  ‘cage’  are  clearly  indicated.  Here  the  output  waveguide  is  shown  on  the  extreme 
right,  its  diameter  as  usual  being  the  same  as  that  of  the  anode  cage. 
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The  septum,  which  is  in  the  form  of  a  circular  disc  with  a  rectangular  cut-out  is  shown  in 
Fig.  18b.  The  thickness  of  the  disc  and  the  height  of  its  rectangular  cut-out  are  both 
fixed,  but  the  width  of  the  cut-out  can  vary  between  10  and  50  mm.  The  three  different 
positions  of  the  disc  relative  to  the  flange  are:  0.0  mm  (right),  40.6  mm  (middle)  and  89.5 
mm  (left).  Thus  we  are  going  to  consider  three  different  positions  of  the  septum  and  four 
different  widths  of  the  rectangular  cut-out,  i.e.  3x4—  12  different  situations,  as  shown 
in  Table  1.5 

1.4.1c  Simulation  results 

An  overview  of  our  simulation  results  is  provided  in  Table  1.5  where  the  peak  output 
power  and  the  corresponding  frequency  of  oscillations  are  shown  as  functions  of  the 
position  of  the  septum  and  the  corresponding  width  of  the  rectangular  cut-out.  In  all 
cases  the  input  power  has  been  kept  constant  at  26  GW,  which  corresponds  to  an  input 
impedance  of  16.7  Q. 
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Table  1.5  Peak  output  power  and frequency  as  functions  of  position  of  the  circular  septum  and  the  width 
of  its  rectangular  act-out. 


Input 

Voltage 

Input 

Current 

Input 

Power 

Peak  Power 

Output 

Frequency 

50mm  Slot  Left 

650  kV 

39  kA 

26  GW 

350  MW 

2.074  GHz 

50mm  Slot  Right 

650  kV 

39  kA 

26  GW 

280  MW 

2.171  GHz 

50mm  Slot  Middle 

640  kV 

40  kA 

26  GW 

175  MW 

2.106  GHz 

30mm  Slot  Left 

640  kV 

40  kA 

26  GW 

200  MW 

2.074  GHz 

30mm  Slot  Right 

650  kV 

39  kA 

26  GW 

150  MW 

2.161  GHz 

30mm  Slot  Middle 

650  kV 

39  kA 

26  GW 

140  MW 

2.051  GHz 

20mm  Slot  Left 

640  kV 

40  kA 

26  GW 

140  MW 

2.081  GHz 

20mm  Slot  Right 

650  kV 

39  kA 

26  GW 

62  MW 

2.174  GHz 

20mm  Slot  Middle 

650  kV 

39  kA 

26  GW 

80  MW 

2.057  GHz 

10mm  Slot  Left 

650  kV 

39  kA 

26  GW 

60  MW 

2.078  GHz 

10mm  Slot  Right 

660  kV 

38  kA 

26  GW 

30  MW 

2.168  GHz 

10mm  Slot  Middle 

650  kV 

39  kA 

26  GW 

58  MW 

2.042  GHz 

1.4.1  cl  Influence  on  the  output  power 

We  find  by  varying  the  width  of  the  rectangular  cut-out,  Fig.  18b,  that  the  best  results  are 
obtained  for  the  biggest  size  of  the  slot  (50  mm  width),  when  the  peak  output  power  is 
the  highest,  the  opposite  being  true  for  the  smallest  size  of  the  slot  (10  mm  width).  Also 
we  can  see  from  Table  1.5  that  placing  the  disc  on  the  left  (i.e.  89.5  mm  away  from  the 
flange)  again  yields  the  best  results;  the  peak  output  power  drops  down  substantially 
when  the  disc  is  moved  near  the  flange  and  it  drops  even  further  when  the  disc  is  placed 
in  its  ‘middle’  position  (40.6  mm  away  from  the  flange),  the  three  different  positions 
being  shown  in  Fig.  1 8a. 
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The  first  of  the  two  effects  is  easy  to  understand,  since  a  small  slot  prevents  the  flow  of 
r.f.  power  towards  the  output,  the  disc  beginning  to  act  as  a  metal  reflector.  The  second 
effect  is  more  difficult  to  unravel  —  it  is  probably  closely  related  to  the  mode  pattern  as 
it  gradually  develops  at  the  build-up  of  oscillations;  one  could  also  think  about  it  as  a 
form  of  impedance  mismatch. 


1.4.1c2  Influence  on  the  frequency  of  oscillations 

Figs.  19a  and  b  respectively  show  the  E-field  and  its  FFT  obtained  in  the  usual  place  i.e. 
near  the  end  of  the  output  waveguide.  Here  we  have  chosen  the  optimum  conditions  of  a 
disc  with  a  50  mm  wide  slot  and  in  its  left  position,  i.e.  89.5  mm  away  from  the  flange, 
Table  1.5. 

At  this  point  it  is  interesting  to  compare  the  different  geometries  of  a  septum  and  the 
corresponding  influence  they  have  on  the  performance  of  a  vircator.  If  we  look  at  Figs. 
13a,  b  which  refer  to  the  rectangular  septum  and  Figs.  16a,  b  which  represent  the  circular 
septum,  then  it  is  clear  that  in  both  cases  the  output  is  much  steadier,  the  corresponding 
frequency  spectrum  being  much  cleaner  than  that  obtained  for  a  vircator  without  a 
septum  of  any  form,  Figs.  3a,  b.  However  the  effect  that  a  septum  has  on  the  resonant 
frequency  seems  to  depend  on  its  geometry.  In  the  case  of  a  rectangular  septum,  the 
resonant  frequency  is  shifted  down  from  2  GHz  to  1.82  GHz;  at  the  same  time  its  value 
strongly  depends  on  the  position  of  the  septum  along  the  axis,  Table  1.4.  In  the  case  of  a 
circular  septum  the  resonant  frequency  is  moved  upwards  from  2  to  2.07  GHz,  its  value 
being  only  slightly  affected  by  the  axial  position  of  the  septum  or  the  size  of  its 
rectangular  slot,  the  total  variation  being  less  than  6  %,  as  can  be  seen  in  Table  1.5. 
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Finally  Figs.  20  and  21  show  the  actual  r.f.  output  power  respectively  for  a  circular  and  a 
rectangular  septum.  We  can  now  see  quite  clearly  that,  although  in  the  presence  of  a 
circular  septum  the  vircator  oscillates  more  smoothly  for  some  60  ns,  its  peak  output 
power  is  only  350  MW;  on  the  other  hand  in  the  presence  of  a  rectangular  septum  the 
peak  output  power  reaches  600  MW  but  only  for  some  20  ns,  the  output  being  somewhat 
less  stable. 

Thus  with  the  aid  of  the  cavity-like  co-axial  geometry  of  the  vircator  and  by  introducing  a 
septum  within  the  anode  mesh,  one  can  achieve  a  reasonably  well  defined  resonant 
frequency  at  the  output  end  of  the  waveguide,  as  shown  in  Figs.  13a,  b  and  16a,  b.  It 
would  appear  that  septa  have  the  effect  of  absorbing  stray  electrons  that  hinder  clean 
virtual  cathode  oscillations.  This  suggestion  is  strongly  supported  by  our  detailed 
investigations  of  electron  trajectories  described  in  Section  1.4.2d. 


1.4.2  Final  model  of  TTU  vircator 
1.4.2a  Introduction 

As  we  pointed  out  before  the  principle  of  generating  microwave  radiation  using  virtual 
cathode  oscillations  is  conceptually  quite  simple  [1,2,4-6,9-15],  however  the  actual 
process  is  highly  non-linear;  it  is  therefore  singularly  difficult  to  treat  it  analytically.  It 
seems  that  computer  modeling  presents  the  only  viable  method  for  the  study  and 
understanding  of  the  vircator  interaction  process.  We  have  therefore  decided  to 
concentrate  in  this  section  on  modelling  in  some  detail  the  cathode/anode  area  where  the 
main  interaction  process  is  taking  place.  We  are  fully  aware  of  the  fact  that  in  the  past 
[1,4,9-14]  some  efforts  have  been  made  to  simplify  matters  by  assuming  that  in  principle 
there  are  two  types  of  interaction  going  on  at  the  same  time.  An  oscillation  due  to 
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forward  and  backward  movement  of  electrons  between  real  and  virtual  cathodes,  the 
semi-permeable  anode  being  situated  in-between  and  a  sort  of  ‘plasma-frequency’ 
oscillation  due  to  the  movement  of  the  virtual  cathode  itself,  mostly  in  the  form  of  space- 
charge  density  fluctuation  around  its  theoretical  saturation  point.  The  two  corresponding 
frequencies,  usually  referred  to  as  freflex  and  fp  ,  invariably  differ  by  a  large  factor  and  any 

attempts  to  use  either  as  an  expected  vircator  frequency  of  oscillations  have  not  been 
very  successful.  The  best  one  could  say  is  that  the  actual  frequency  of  oscillations  occurs 
in  the  following  range  [1,4,10] 

2-2 /p  <fosc  <2.7 4 

We  will  see  more  clearly  later  on  some  reasons  for  the  above  discrepancy. 

A  close  observation  of  the  interaction  process  as  displayed  by  the  computer  model 
revealed  that,  as  perhaps  we  should  have  expected,  the  movement  of  electrons  is  much 
more  complex  than  it  was  generally  assumed.  Certainly  in  the  case  of  a  coaxial  vircator 
there  is  a  strong  whirl-like  motion  around  the  metal  ring  situated  at  the  end  of  the  anode 
cage,  Figs.  27a-f.  This  motion  clearly  affects  the  frequency  at  which  the  virtual  cathode 
oscillates  around  its  saturation  point.  In  addition  we  have  a  marked  axial  component  of 
current,  both  inside  and  outside  of  the  anode  cage.  In  fact  the  electrons  forming  the 
space  charge  appear  to  be  divided  into  three  groups: 

•  those  intermittendy  forming  a  virtual  cathode  on  the  inside  of  the  anode, 

•  those  in  axial  circular  motion  hugging  the  anode  or  whirling  around  the  support 
ring  and 


(1.1) 
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•  the  ‘rogue’  or  stray  electrons  which  move  towards  the  output  waveguide  and  thus 
remove  energy  from  the  interaction  mechanism  and  lower  electronic  efficiency  of 
the  power  conversion  process. 

Having  obtained  some  insight  into  the  complexity  of  the  interaction  process  we  were 
then  able  to  improve  the  efficiency  of  the  device  by  trying  to  get  rid  of  the  rogue 
electrons.  As  explained  in  Section  1.3.3d  this  was  achieved  by  altering  the  design  of  the 
anode  and  introducing  a  carefully  shaped  septum  placed  inside  the  anode,  Figs.  28a  and 
b.  The  septum  succeeds  in  trapping  the  stray  electrons  which  intern  increases  the 
microwave  power  output  by  some  52  %. 


1.4.2b  Input  and  output  pulses 

As  a  first  step  in  our  simulation  process  we  have  decided  to  repeat  the  results  shown  in 
Figs.  22  and  23.  Introducing  absorbing  boundaries  at  both  ends  of  our  model  and 
applying  500  kV  pulse  of  70  ns  duration,  we  have  obtained  simulated  input  voltage  and 
current  as  functions  of  time,  as  shown  in  Fig.  24.  The  peak  values  of  730  kV  and  38.4  kA 
are  within  1  -  5  %  of  the  measured  values  shown  in  Fig.  22.  The  contour  of  the 
measured  pulse  is  much  more  rugged  than  that  obtain  by  simulation,  but  this  tends  to  be 
true  of  many  simulation  processes,  the  rugged  shape  being  influenced  by  local 
construction  details  not  included  in  the  model  (e.g.  the  thin  rods  supporting  the  anode). 
Again  dividing  peak  voltage  by  peak  current  we  obtain  19.0  Q,  which  is  within  5  %  of  the 
measured  value. 

Let  us  now  consider  the  simulated  input  and  microwave  output  power,  as  shown  in  Fig. 
25;  we  find  that  their  respective  peak  values  are  25.8  GW  and  600  MW,  Poynting’s 
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theorem  being  used  to  obtain  the  simulated  output  microwave  power.  The  computed 
peak  values  are  again  within  5-8  %  of  the  measured  values  shown  in  Fig.  23.  Both  sets  of 
results  show  two  distinct  peaks  in  the  microwave  output  power;  this  could  well  be  due  to 
some  impedance  mismatch  in  the  dynamics  of  virtual  cathode  oscillations.  Also  the  range 
of  small  fluctuations  in  the  computed  output  is  typical  of  mode  overlap  in  the  heavily 
overmoded  output  waveguide.  We  find  from  Fig.  23  by  dividing  peak  output  by  peak 
input  power  that  the  electronic  efficiency  of  the  tube  is  of  the  order  of  2.3  %,  which  is 
usual  for  vircators.  The  corresponding  simulation  value  obtained  from  Fig.  25  is  also  2.3 
%  in  close  agreement  with  the  measured  value.  This  further  supports  our  confidence  in 
the  adequacy  of  our  modelling  process. 

Finally  we  note  that  there  is  a  natural  time  delay  between  the  input  pulse  measured  near 
the  cathode/anode  space  and  the  output  power  measured  near  the  end  of  the  waveguide, 
Figs.  23  and  25.  In  the  computer  model  shown  in  Fig.  25,  the  delay  is  approx.  13  ns  and 
it  can  be  accounted  for  by  the  time  it  takes  the  pulse  to  travel  down  the  system. 
(Analytical  expressions  for  the  time  of  travel  cannot  be  easily  obtained  because  the 
output  waveguide  is  heavily  overmoded).  The  measured  delay,  Fig.  23,  is  somewhat 
longer  and  of  the  order  of  20  ns,  probably  due  to  a  slower  build  up  of  oscillations  than 
predicted  by  the  model. 

1.4.2c  Frequency  spectrum 

The  simulated  frequency  spectrum  is  shown  in  Fig.  26.  It  was  obtained  by  applying  FFT 
to  the  output  from  an  E-field  probe  situated  near  the  window  of  the  output  waveguide. 
We  find  that  the  main  resonant  frequency  occurs  at  2.108  GHz,  which  compares 
favourably  with  the  measured  value  of  2.00  GHz  (this  is  a  different  resonant  frequency 
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than  that  quoted  in  Fig.  10b,  the  difference  being  due  to  more  recent  TTU  structural 
adjustments  concerning  the  cathode/ anode  gap,  which  we  included  in  our  model).  The 
frequency  spectrum  is  still  fairly  broad-band,  in  spite  of  the  presence  of  a  rectangular 
septum,  Figs.  6  and  7a,  b  -  this  situation  is  typical  of  vircators.  Without  a  septum  the 
frequency  spectrum  is  quite  chaotic,  as  shown  in  Fig.  3b. 


1.4.2d  Electron  trajectories 

In  order  to  consider  the  interaction  process  in  detail  it  has  been  decided  to  create  a 
movie  which  would  simultaneously  trace  the  movement  of  some  10,  000  electrons  as 
they  are  emitted  by  the  cathode,  approach  or  pass  the  foil  anode  and  then  form  a  virtual 
cathode.  The  MAGIC  code  we  have  used  is  ideally  suited  for  this  kind  of  investigation. 
We  then  discovered  that  instead  of  reflexing  radially  between  real  and  virtual  cathodes,  as 
postulated  by  some,  the  movement  of  the  electrons  is  much  more  complicated.  First  of 
all  there  is  a  short  period  lasting  some  6  ns,  when  the  electron  cloud  is  being  set  up,  the 
electron  cloud  remaining  relatively  uniform,  as  shown  in  Fig.  27a.  We  then  find  that  the 
electrons  split  into  three  distinct  groups.  The  first  group  does  indeed  pass  through  the 
semipermeable  anode  and  forms  a  virtual  cathode  which  then  saturates,  sending  electrons 
back  towards  the  cathode,  the  repetition  of  this  process  giving  rise  to  the  usual  form  of 
relaxation  oscillations.  However  not  all  electrons  take  part  in  this  process.  A  second 
group  of  electrons  fails  to  contribute  to  the  virtual  cathode  and  either  travels  in  the  axial 
direction  close  to  the  anode,  Fig.  27b  or  even  wraps  itself  in  tight  circles  round  the  anode 
ring,  Fig.  27c.  Finally  there  is  a  third  group  of  electrons,  called  ‘stray’  or  ‘rogue’  electrons, 
which  detach  themselves  from  the  cloud  and  travel  in  bunches  down  the  system,  finally 
landing  either  on  the  anode  near  the  flange  or  even  on  the  walls  of  the  waveguide,  Figs. 
27d.  All  this  can  be  seen  very  clearly  in  the  ‘movie’  obtained  by  observing  Figs.  27b-f  in 
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quick  succession  on  the  screen.  The  stray  electrons  detract  from  the  proper  interaction 
process  and  should  be  eliminated  if  at  all  possible.  Figs.  27b-f  show  a  complete  cycle  of 
the  interaction  process  which  lasts  for  T  =  0.52  ns,  which  corresponds  to  the  f  —  1/T  = 
1.923  GHz  and  is  within  9  %  of  the  oscillation  frequency  £=2.108  GHz,  Fig.  26.  The 
exact  value  should  have  been  T  =  0.47  ns,  but  it  was  difficult  to  freeze  the  ‘frame’  which 
would  correspond  more  closely  to  the  end  of  one  cycle  of  oscillations. 


1.4.2e  Suggestion  for  an  improved  design 

Having  unravelled  with  the  help  of  Figs.  27a-f  some  of  the  intricacies  of  the  cylindrical 
vircator  interaction  process,  we  have  decided  to  investigate  the  possibility  of  improving 
the  performance  of  the  tube.  An  obvious  target  was  to  get  rid  of  the  so-called  stray 
electrons  mentioned  in  the  previous  section.  About  eight  or  ten  septa  of  varying 
geometry  have  been  investigated,  until  finally  the  design  shown  in  Figs.  28a  and  b  proved 
to  be  the  best.  In  addition  it  was  necessary  to  replace  the  semipermeable  mesh  by  solid 
metal  over  one  half  of  its  length  near  the  flange,  as  shown  in  Fig.  28a.  All  these  changes 
required  many  additional  computations,  each  run  lasting  for  some  12  —  16  horns. 

The  input  and  output  power  pulses  are  shown  in  Fig.  29.  The  model,  shows  that  for  the 
same  peak  input  power  of  25.9  GW,  the  peak  output  power  has  grown  from  600  to  913 
MW,  an  increase  of  52  %.  At  the  same  time  the  electronic  efficiency  calculated  as  the 
ratio  of  peak  output  to  peak  input  power  has  grown  from  2.3  %  to  3.5  %,  which  is  quite 
a  welcome  improvement. 

Finally  the  FFT  of  the  output  from  an  E-field  probe  situated  near  the  end  of  the 
waveguide  is  shown  in  Fig.  30.  Comparing  Figs.  30  and  26  we  find  that  the  main 
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resonance  frequency  has  shifted  upwards  from  2.108  to  2.153  GHz,  the  whole  spectrum 
being  similar  to  that  shown  in  Fig.  26.  The  shift  in  frequency  was  to  be  expected  since  it 
has  been  noted  before  that  the  main  resonant  frequency  is  very  sensitive  to  the  thickness 
of  the  anode  ring  and  the  length  of  the  permeable  section  of  the  anode.  As  far  as  the 
‘cleaning  up’  of  the  spectrum  is  concerned  we  noted  on  several  occasions  that  septa  have 
that  effect,  presumably  by  acting  somewhat  like  a  quasi-cavity.  Here  one  has  to  strike 
balance  between  this  effect  and  a  drop  in  power  caused  by  the  septum  acting  as  an 
obstruction  to  the  flow  of  power. 
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Conclusions  -  part  1. 

Two  main  conclusions  can  be  drawn  from  the  above  investigations.  In  the  first  place  we 
have  established  that  computer  modelling  of  vircators  has  reached  a  degree  of 
sophistication  consistent  with  its  use  as  a  design  tool.  Secondly  we  were  able  to  show  that 
the  interaction  process  between  the  beam  and  the  e.m.  field  is  not  only  highly  nonlinear 
but  also  very  complex,  so  that  an  analytical  approach  would  be  fraught  with  difficulties, 
computer  modelling  being  the  only  practical  solution.  Finally  having  investigated  the 
interaction  process  in  some  detail  we  have  been  able  to  suggest  changes  in  design  which 
should  increase  the  microwave  power  output  by  some  52%. 
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Figures  -  part  1 


Figure  1  Schematic  diagram  of  a  plane  vircator. 

Figure  2  Cross-section  of  an  early  model  of  the  TTU  vircator. 

Figure  3a  Ep  field  intensity  as  recorded  near  the  output  window. 

Figure  3b  Frequency  spectrum  of  the  microwave  output  (FFT  of  Fig.  3a)  —  main 
resonant  peak  at  2.73  GHz. 

Figure  4  Snapshot  of  Ep  contours  near  the  output  window. 

Figure  5  Snapshot  of  Bp  contours  near  the  output  window. 

Figure  6  Gun  Section  of  the  TTU  vircator. 

Figure  7a  Axial  cross-section  of  the  improved  TTU  vircator  model. 

Figure  7b  Transverse  view  of  the  rectangular  septum  situated  in  plane  A-A’  of  Fig.  6a. 
Figure  8a  Rectangular  septum  and  part  of  the  anode  in  Cartesian  co-ordinates. 

Figure  8b  Rectangular  septum  and  part  of  the  anode  in  polar  co-ordinates. 

Figure  9  Various  ways  of  applying  the  anode  voltage. 

Figure  10a  Typical  vircator  output  -  TTU  E-field  probe  measurements. 

Figure  10b  Frequency  spectrum  —  fast  Foruier  transform  (FFT)  of  Fig.  10a  —  main 
resonant  peak  at  1.834  GHz. 

Figure  1 1  Simulated  input  power  pulse. 

Figure  12  Simulated  microwave  output  power  pulse. 

Figure  13a  Simulated  E-probe  output. 

Figure  13b  FFT  of  Fig.  13a  -  peak  resonant  frequency  at  1.816  GHz. 

Figure  14  Resonant  frequency  as  a  function  of  cathode/anode  gap. 

Figure  15  Peak  resonant  frequency  as  a  function  of  axial  position  of  the  rectangular 
septum,  Fig.  8b. 

Figure  1 6  Enlarged  smooth  edge  septum. 

Figure  17  Frequency  spectrum  with  enlarged,  smooth-edge  septum  in  place. 

Figure  18a  Axial  Cross  Section  of  the  interaction  region. 

Figure  18b  Septum  in  the  form  of  a  circular  disc  with  a  rectangular  cut-out. 

Figure  19a  E-field  output  in  the  presence  of  a  circular  septum  (‘lcft’position,  slot  50mm 
wide) 

Figure  19b  FFT  plot  of  Fig.  19a;  main  resonance  now  at  2.074  GHz. 
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Figure  20  Output  power  pulse  in  the  presence  of  a  circular  septum  (‘left’  position,  slot 
50mm  wide) 

Figure  21  Output  power  pulse  in  the  presence  of  a  rectangular  septum. 

Figure  22  Typical  TTU  input  voltage  and  current  pulses. 

Figure  23  Typical  TTU  input  power  and  microwave  output  power  pulses. 

Figure  24  Simulation  input  voltage  and  current  pulses. 

Figure  25  Simulation  input  and  microwave  output  power  pulses. 

Figure  26  Simulated  frequency  spectrum  —  main  resonant  frequency  at  2.108  GHz. 

Figure  27a  Electron  cloud  prior  to  the  start  of  oscillations,  t  =  6.60  ns. 

Figure  27b  Electron  cloud  at  t  =  20.00  ns;  arrow  indicating  the  electrons  lingering  near 
the  anode  foil  and  drifting  towards  the  flange. 

Figure  27c  Electron  cloud  at  t  =  20.12  ns;  arrow  indicating  electgrons  curling  around  the 
anode  supporting  ring. 

Figure  27d  Electron  cloud  at  t  =  20.24  ns;  arrow  indicating  electrons  moving  in  bunches 
towards  the  waveguide. 

Figure  27e  Electron  cloud  at  t  =  20.36  ns. 

Figure  27f  Electron  cloud  at  the  end  of  one  cycle  of  oscillations,  t  —  20.52  ns. 

Figure  28a  Cross-section  of  the  vircator-improved  design. 

Figure  28b  New  septum  design;  plane  B-B’  of  Fig.  28a. 

Figure  29  Simulation  input  and  microwave  output  power  pulses  -  improved  design. 
Figure  30  Frequency  spectrum  —  improved  design;  main  resonance  at  2.153  GHz. 
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Figure  1  Schematic  diagram  of  a  plane  vircator. 
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Figure  2  Cross-section  of  an  early  model  of  the  TTU  vircator. 
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Erho  at  PROBE. Z1 40 .R9 


Figure  3a  Ep  field  intensity  as  recorded  near  the  output  window. 


Magnitude  Of  FFT  Of  Erho  at  PROBE. Z140.R9 


Figure  3  b  Frequency  spectrum  of  the  microwave  output  (FFT  of  Fig.  3  a)  -  main  resonant  peak  at  2.73 


GHZ. 
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Time  7.2B2  ns:  Of  Eiho  (V/m) 


Z  (H)  (E_Z) 


(a)  early  stage  of  oscillations,  t— 7.282  ns. 


(b)  advanced  stage  of  oscillations,  t—  1 1.029  ns. 

Figure  4  Snapshot  of  Ep  contours  near  the  output  window. 
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Time  8.093  ns:  Of  Biho  (tesla) 
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(a)  early  stage  of  oscillations,  t-8.093  ns 


Time  11.029  ns:  Of  Biho  (tesla) 
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Figure  5  Snapshot  of  Bp  contours  near  the  output  window. 


44 


1.17  1.18  1.19  1.2  1.21  1.22  1.23  1.24  1.25  1.26 
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Figure  10a  Typical  vircator  output  —  TTU  E-field  probe  measurements. 


Scope  Output 


Figure  10b  Frequency  spectrum  -  fast  Fourier  transform  (FFT)  of  Fig.  10a  -  main  resonant  peak  at 
1.834  GHz- 
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Frequency  (GHz) 

Figure  13b  FFT  of  Fig.  13a—  peak  resonant frequency  at  1.816  GH 
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Number  of  Protruding  Emitter  Cell 


Peak  Resonant  Frequency  Response  vs  Number  of  Protruding  Emitter  Cell 


Peak  Resonant  Frequency  (GHz) 


Figure  14  Resonant frequency  as  a  function  of  cathode /  anode  gap. 
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Figure  16  Enlarged  smooth  edge  septum. 


8a  Axial  Cross  Section  of  the  interaction  region. 
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Figure  19a  E-field  output  in  the  presence  of  a  circular  septum  ('left’ position,  slot  50mm  wide) 


FIELD  Magnitude  Of  FFT  Of  Erho  at  RBB . P90 .Z1500 


Figure  19b  FFT plot  of  Fig.  19a;  main  resonance  now  at  2.074  GH £ 
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FIELD  POWER  S.DA  at  OUT. RIGHT. IN 


Figure  20  Output  power  pulse  in  the  presence  of  a  circular  septum  (‘left’ position,  slot  50mm  wide) 


FIELD_POWER  S.DA  at  OUT. RIGHT. IN 


Figure  21  Output  power pulse  in  the  presence  of  a  rectangular  septum. 
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Figure  23  Typical  TTU  input  power  and  microwave  output  power pulses. 


Figure  25  Simulation  input  and  microwave  output  power  pulses. 
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Figure  27 d  Electron  cloud  at  t  =  20.24  ns;  arrow  indicating  electrons  moving  in  bunches  towards  the 
waveguide. 


1  (in) 

Figure  27 f  Electron  cloud  at  the  end  of  one  cycle  of  oscillations,  t  —  20.52  ns. 
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Figure  28a  Cross-section  of  the  vircator-improved  design. 


Figure  28b  New  septum  design; plane  B-B ’of  Fig.  28a. 


Figure  30  Frequency  spectrum  -  improved  design;  main  resonance  at  2. 153  GHp. 


Part  2. 


2.1.  Introduction 

In  principle  the  electron  beam  generating  element  in  an  axial  vircator,  such  as  that  shown 
in  Fig.  1,  can  be  represented  by  an  inverted  cylindrical  diode  with  the  cathode  on  the 
outside  and  the  anode  on  the  inside  [1].  In  the  case  of  a  vircator  it  is  common  to  use  a 
velvet  cathode  which  is  capable  of  delivering  very  high-density  electron  emission  for  a 
period  of  ~  100  ns.  Since  the  anode  is  semi-permeable,  the  electron  beam  passes  through 
it  virtually  unhindered  and  enters  the  space  inside  the  anode,  where  a  virtual  cathode  can 
be  formed,  Fig.  2.  It  is  this  part  of  the  device  which  we  are  now  going  to  consider. 

In  order  to  simplify  the  algebra  of  the  problem  we  choose  plane  rather  than  cylindrical 
geometry  and  use  a  short-circuited  plane  diode  as  a  model  for  the  space  within  the 
vircator’s  cylindrical  anode,  see  II  in  Fig.  3.  Here  we  assume  that  a  single-velocity 
electron  stream  is  generated  in  I  and  enters  the  inter-electrode  space  II,  the  entrance 
velocity  of  the  electrons  being  governed  by  the  anode/ cathode  potential  difference  Va  . 
Since  the  diode  II  is  short-circuited,  the  exit  electrode  (collector)  is  kept  at  the  same 
potential  as  the  entrance  electrode  (anode).  Also  it  is  assumed  that  most  of  the  time  the 
anode  collects  no  current,  its  sole  purpose  being  to  establish  an  idealised  plane  at 
potential  Va  relative  to  the  cathode.  In  the  actual  vircator  there  is  no  collector  and  the 
electrons  gradually  lose  their  radial  velocity  as  they  approach  the  axis,  as  shown  in  Fig.  2, 
but  we  would  like  to  suggest  that  the  addition  of  the  collector,  which  is  necessary  to 
complete  the  circuit,  is  not  likely  to  affect  the  basis  of  our  argument. 
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Those  of  us  who  stdll  recall  vacuum  tubes  may  remember  that  a  system  shown  in  Fig.  3 
was  thoroughly  analysed  by  no  lesser  a  man  than  Shockley  [2].  However,  he  and  others 
did  it  in  the  absence  of  a  magnetic  field  and  our  task  will  be  to  see  what  happens  when 
we  introduce  an  external  magnetic  field  parallel  to  the  electrodes,  as  indicated  in  Fig.  3. 
Here,  we  have  limited  the  presence  of  the  magnetic  field  to  region  II;  its  presence  in 
region  I  can  be  safely  ignored,  since  no  virtual  cathode  would  be  formed  there,  region  I 
merely  acting  as  a  source  of  current  for  region  II. 
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2.2.  Magnetic  diode  model 

It  is  generally  recognised  that  the  presence  of  a  magnetic  field  substantially  increases  the 
power-handling  capacity  of  a  microwave  tube;  this  feature  greatly  contributed  to  the 
success  of  the  so-called  M-type  devices  [3].  A  simple  physical  explanation  of  the 
phenomenon  can  be  gained  with  the  help  of  Fig.  4,  where  we  have  a  graphic 
representation  of  the  Poisson  equation  in  the  absence  and  in  the  presence  of  a  magnetic 
field.  In  the  absence  of  the  field  we  have  a  ‘simple’  form  of  the  Poisson  equation  given 
by  (2.1)  and  shown  in  Fig.  4,  where  curve  (1)  is  a  hyperbola  asymptotic  to  horizontal  and 
vertical  axes. 

^l  =  -n  =  -  I  =  __L 

dx2  x" 

Here  we  have  used  the  ‘bar’  notation  as  explained  in  Appendix  1  in  order  to  simplify  the 
algebraic  form  of  the  equations.  The  main  advantage  of  the  bar  notation  is  that  the 
distance  between  the  electrodes  x  =  x/d  ,  and  the  electron  entrance  velocity  v  =  v/ Vj  are 
both  equal  to  unity. 


2  _ 

In  the  absence  of  the  magnetic  field  x'  =<f>  ,  which  leads  direcdy  to  (2.1).  The 

introduction  of  a  magnetic  field  alters  the  situation  substantially  since  now  y'  ^  0  and 

_  2  2 

hence  the  new  energy  balance  is  expressed  by  (j)  =x'  +y*  (see  (2.13)).  Substituting  for 
X-  in  the  Poisson  equation  we  now  obtain 
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The  corresponding  curve  (2)  in  Fig.  4  is  asymptotic  not  to  the  vertical  axis  ^  =  0  ,  but  to 
a  vertical  line  (f>  =  cy20X2  >0  and  thus  lies  everywhere  above  curve  (1).  The  physical 
meaning  of  the  difference  between  curves  (1)  and  (2)  is  quite  significant.  Let  us  consider 
the  conditions  at  a  fixed  potential  (j)  and  a  fixed  injection  current  density  Jc .  Fig.  4 
indicates  that  in  the  presence  of  the  magnetic  field  the  system  can  absorb  a  much  higher 
space-charge  density  p  ,  curve  (2)  lying  everywhere  above  curve  (1),  the  difference 

growing  as  we  approach  (j)  —  ft>20x2 .  Simplistically  one  could  say  that  the  electrons  now 
have  more  ‘room’  for  manoeuvre,  being  able  to  move  sideways  as  well  as  forward. 

Algebraically  the  difference  relies  on  the  fact  that  only  x‘  enters  in  the  calculation  of  the 
injection  current  density  Jc ,  whereas  the  energy  balance  expressed  by  <f>  incorporates 

both  X'  and  y*. 
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2.3  Solution  in  terms  of  time 


It  is  well  known  that  in  the  analysis  of  steady-state  electron  beams  it  is  algebraically  much 
simpler  to  use  time  t  rather  than  distance  x  as  an  independent  variable,  the  two 
variables  being  interchangeable,  x  =  x(t) .  We  propose  therefore  to  carry  out  our  analysis 
in  terms  of  t  . 


2.3.1  Basic  equations 

Using  our  bar  notation  the  usual  equations  of  motion  for  a  single  electron  in  the 
presence  of  a  magnetic  field  are  given  by: 


x" 


-Ie- 

2 


®coy‘ 


(2.3) 


y"  =  ®cox' 


(2.4) 


where  E  =  — r  ,®co  =eB/m0  and  O)c0  is  defined  in  Appendix  1.  Similarly  in  a  steady  state 
dx 

the  current  density  is  defined  by  writing: 


Jc 


E- 


(2.5) 


Differentiating  (2.3)  once  with  respect  to  time,  substituting  from  (2.4)  and  then  again 
from  (2.5),  we  obtain 
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x"  =--£•■ 
2 


•®coy' 


1  -7  — 2  — . 

—  Jo 


or 


x-  +  (Wc20x‘ 


=  const 


(2.6) 


Jc  being  constant  in  ID.  We  can  easily  recognise  (2.6)  as  the  DC  part  of  the  Llewellyn- 
Peterson  equation  when  a  magnetic  field  is  present;  the  effect  of  the  space  charge  is 

in  (2.6)  by  the  current  density  term  -  ^  Jc  >  0  ,  electrons  carrying  negative 
charge.  The  solution  of  (2.6)  is  quite  straightforward  bearing  in  mind  that  in  our  notation 
we  have  at  the  entrance  electrode:  Xj  =  0 ,  xi  =  1  and  yj  =  0  . 


^  =  ^i-(1-c)  +  1i-s--^(0-s)  =  --^(1-c)  +  ^s  +  -^(0-s)  (2.7) 


1  2i 


®c0  ®c0  1(0 


cO  2fi)c0 


®c0  9©c0 


X-  =-^i-S  +  C - ir-(l-c) 

<yco  2<yco 


Ei  2 i  t  \ 

—  S  +  C  +  — ^(l-c) 


x”  =x;'c-fflr0s - - — s 

2®c0 


where 


2®c0  9 0cO 


E,  _  2 1 

= - -  C  —  ®cnS  + - S 

2  9*>c0 


(2.8) 


(2.9) 


t  =  t-tj ,  s  =  sin«yc0r  ,  c  ~  cos 0)c0t  ,  6  =  <»c0r 


(2.10) 


72 


and 


from  (2.3) 


(2.11) 


from  (A  1.2) 


(2.12) 


Once  E,  is  known  for  any  give  l ,  the  problem  of  potential  and  velocity  distribution  is 

solved  since  x,  x'  can  be  obtained  from  (2.7),  (2.8)  and  y,  y*  from  integrals  of  (2.4). 
This  is  a  much  quicker  procedure,  especially  in  the  presence  of  an  initial  velocity,  than 
any  attempt  to  solve  the  Poisson  equation  (2.2)  directly,  as  was  common  in  the  past  for  a 
zero  initial  velocity  [4,5]. 

At  this  point  it  is  convenient  to  introduce  the  usual  energy  balance  equation.  Multiplying 
(2.3)  by  x‘ ,  (2.4)  by  y' ,  adding  and  then  integrating  with  respect  to  t  we  obtain 

</>  =x'  +y‘  (2-13) 

Here  we  have  assumed,  bearing  in  mind  Fig.  3,  that  all  electrons  start  with  zero  velocity 
at  a  cathode  which  is  at  zero  potential,  so  that  <f>c  =  X'c  =  y|.  =  0 .  Thus,  at  the  entrance 

electrode  we  have,  by  definition  (see  Appendix  1):  Xj  =  0 ,  $  =  Xj  =  1 ,  yj  =  0  .  Since  by 
integrating  (2.4)  once  we  obtain 


73 


(2.14) 


y=®cox 


the  substitution  of  (2.14)  in  (2.13)  yields 


<f  +<yc20x2 


(2.15) 


For  a  short-circuited  diode  we  have  at  the  exit  electrode  x2  =  <j)2  =  1  and  the  energy 
balance  (2.15)  reduces  to 


x-  =l-tyc0 


(2.16) 


Soon  we  will  use  (2.16)  as  one  of  the  boundary  conditions. 


2.3.2  Derivation  of  E,  (?) 

The  next  step  in  our  investigations  is  to  obtain  a  relationship  between  the  injection 
current  density  l  and  the  electric  field  at  the  entrance  electrode  Ej  .  In  the  absence  of  a 

magnetic  field  this  relationship  can  be  expressed  directly  in  the  form  l  =  l (Ej )  by  writing 
(see  (A2.12)  in  Appendix  2) 


y =E, 


(2.17) 
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This  is  no  longer  possible  in  the  presence  of  the  magnetic  field.  The  best  we  can  do  now 
is  to  obtain  the  relevant  curve  in  its  parametric  form  by  writing  *(rd)  ,  where  fd 

is  the  total  transit  time  between  the  electrodes.  By  assuming  yj  =  0 ,  putting  x  =  x2  =  1  in 

(2.7)  and  x'2  =  \Jl  - in  (2-8),  we  obtain  from  the  two  equations: 


2F52  ^  Cd^  +  f)  S<l  +  ()m3  ^d  Sd^ 
2tyc0  coc0  9tyc0 


(2.18) 


x; 


=r- 


co, 


E, 


2 1 


cO 


Sa+cd+— r(l-cd) 
2coc0  9tyc0 


(2.19) 


where 


sd  =  sintyc0rd ,  cd  =  cos«c0rd ,  0d  =  coc0Td 


(2.20) 


Here  as  before  the  subscripts  1  and  2  refer  respectively  to  entrance  and  exit  electrodes. 
Substituting  for  E,  from  (2.19)  we  obtain  from  (2.18),  after  some  heavy  algebra 


(2.21) 


We  can  now  substitute  (2.21)  in  (2.19)  in  order  to  obtain  the  following  expression  for  the 
electric  field  at  the  entrance  electrode: 
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Ej  —  2©c0 


2i 


tan 


(2.22) 


|l-7l-©c0 

~(l  +  Vl”®cO  j 

tan2 

fl  _  _  ) 
7 

2  tan 

n_ 

(\ _ 

_  ®c0^d 
\2 


At  the  cut-off  we  have,  by  definition,  \'2  =  yjl - ®c20  =0,  so  that  a>c0  =1  and  (2.21),  (2.22) 
acquire  a  simpler  form  given  by 


(2.23) 


(2.24) 


An  expression  similar  to  (2.21)  has  been  derived  by  others  [6],  but  their  boundary 
conditions  and  notation  differ  from  ours. 

A  graphical  presentation  of  the  results  is  provided  in  Fig.  5  where  we  have  the  electric 
field  at  the  entrance  electrode  Ej  as  a  function  of  the  injection  current  density  i  for 
different  values  of  the  magnetic  field  B,  expressed  in  terms  of  the  reduced  electron  gyro- 
frequency  <ac0  defined  in  Appendix  1.  Two  points  of  particular  interest  have  been 
respectively  marked  A  and  B  in  Fig.  5  and  they  will  be  considered  separately.  It  should  be 
noted  that  all  curves  start  at  the  origin  in  Fig.  5;  this  confirms  the  obvious  fact  that 
irrespective  of  the  value  of  the  magnetic  field,  the  electric  field  at  the  entrance  electrode 
Ej  =  0  in  the  absence  of  space  charge,  p  =  0  or  z  =  0  ina  single-stream  conducting  diode. 
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Point  A 


At  point  A  the  injection  current  density  l  reaches  its  maximum  value  for  a  given  coc0 . 
The  value  /max  is  obtained  by  differentiating  (2.21)  with  respect  to  rd  and  then  putting 

=  0  .  After  some  algebra  we  obtain 
drd 

(l  +  Vl-^co)(<9d  "Sd)  =  ^cO  (l-Cd)  (2-25) 


which  must  be  satisfied  by  Td  for  any  given  O)c0 .  The  maximum  value  of  the  injection 
current  is  then  given  by 


2m\, 

0d  ~sd  l-cd 


(2.26) 


Here  6d  ,  sd  and  cd  are  defined  in  (2.20). 

We  now  find  from  Fig.  5  that  as  the  magnetic  field  increases,  i.e.  as  a>c0  grows  from  zero 
for  B  =  0  to  G)c0  =  1  for  the  cut-off  value  B  =  BH ,  the  maximum  value  of  the  injection 
current  l  (point  A)  decreases  until  it  suddenly  drops  to  zero  for  B  >  BH  (not  shown).  All 
this  is  in  agreement  with  our  discussion  of  the  general  effect  of  the  magnetic  field 
presented  in  Section  2.2.  As  the  magnetic  field  increases  at  a  fixed  value  of  the  space 
charge  density  p  the  tangential  component  -Jy  grows  at  the  expense  of  the 

perpendicular  component  —  Jx  =  —  Jc  =(4/9)/  ;  thus  for  B  =  0  or  G)c0  —  0  we  have 
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or  l  =  0  does  not  imply  p  =  0  in  the  presence  of  a  double-stream  flow).  Table  2.1  gives 
the  values  of  zmax  and  ,5^  =  (4/9) /max  for  several  values  of  the  magnetic  field  expressed 
by  (Oc0  .  The  curve  joining  all  the  /  =  lmax  or  A  points  is  shown  dotted  in  Fig.  5.  For  the 
sake  of  completeness  the  corresponding  values  of  the  transit  time  rd  max  ,  the  transit  angle 
0msx  =  d)c0fd  max  and  the  electric  field  at  the  entrance  electrode  Ej  have  also  been 
included. 

Typical  potential  (f>  and  velocity  X‘2  distributions  for  different  values  of  <S>c0  are  shown 
in  Figs.  6a  and  b  respectively.  Here  for  typographic  convenience  we  have  called  the  x- 

directed  velocity  u  rather  than  X' . 

Point  B 

Let  us  now  consider  the  second  point  of  interest,  viz  point  B.  It  has  been  shown  by 
others  [7]  that  in  the  absence  of  the  magnetic  field,  COc0  =  0  ,  all  points  lying  on  the  lower 
branch  of  the  curve  between  the  origin  and  point  A  represent  stable  conditions,  whereas 
all  points  between  A  and  B  represent  unstable  solutions,  the  curve  becoming  meaningless 
beyond  B  due  to  the  formation  of  a  virtual  cathode  which  introduces  discontinuity  in  l . 
Physically  when  an  attempt  is  made  to  increase  the  injection  current  l  beyond  its 
maximum  value  tmax ,  the  system  rapidly  collapses:  first  the  potential  j>  and  hence  the 
electron  velocity  precipitously  drops  to  zero,  point  B  in  Fig.  5  and  then  the  position  of 
the  virtual  cathode  is  shifted  by  splitting  l  into  l2  reaching  the  exit  electrode  and 
^  =  l  — 12  returning  to  the  entrance  electrode.  The  best  physical  description  of  the 


process  can  be  found  in  [8].  It  is  now  suggested  that  similar  conditions  should  prevail  in 
the  presence  of  the  magnetic  field,  i.e.  that  the  lower  section  of  each  curve  for 
0  <  i  <  imax  would  represent  stable  solutions,  whereas  the  upper  part  of  each  curve  would 

indicate  unstable  solutions.  Since  the  operation  of  a  vircator  depends  on  the  formation  of 
a  virtual  cathode,  it  is  essential  to  know  what  happens  to  it  in  the  presence  of  a  magnetic 
field,  bearing  in  mind  that  now  <f>  is  no  longer  equal  to  x'2  but  is  given  by  (2.15). 

At  this  point  it  is  essential  to  distinguish  between  two  types  of  instability.  In  Fig.  9  curve 
A'  is  stable,  curve  A  is  conditionally  stable  (l  =  zmax)  and  curve  B  is  inherently  unstable; 
that  means  that  in  the  case  of  curve  B  the  real  part  of  the  corresponding  Llewellyn- 
Peterson  coefficient  (in  this  case  an)  has  a  positive  real  part.  In  this  case  of  curves  A' 
and  A  the  real  part  of  an  is  respectively  negative  and  zero.  When  we  attempt  to 
increase  the  injection  current  l  beyond  its  maximum  value  tmax ,  the  system  jumps  from 
solution  A  ,  to  solution  C  ,  momentarily  going  through  a  potential  distribution 
resembling  B  (the  corresponding  idealized  ‘stead  state’  solution  would  have  l  <  Zmax  ). 

The  solution  C  is  characterised  by  the  appearance  of  a  virtual  cathode  and  the 
corresponding  partial  reflection  of  current.  Although  possible  in  theory,  such  a  situation 
cannot  exist  in  practice,  largely  due  to  inhomogeneities  of  the  system.  Therefore  in 
practice  the  solution  becomes  unstable,  the  system  continuously  relaxing  between 
solutions  with  and  without  a  virtual  cathode,  thus  giving  rise  to  oscillations.  Exactly  the 
same  argument  applies  in  the  presence  of  the  magnetic  field,  as  indicated  in  Figs.  10a  and 
b,  Fig.  8a  merely  indicating  a  general  scheme  of  things.  What  is  important  from  the  point 
of  view  of  a  vircator,  that  as  the  magnetic  field  increases,  the  amplitude  of  the 
corresponding  oscillations  decreases  until  they  are  finally  extinguished  when  the  magnetic 
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field  approaches  its  cut-off  value.  In  the  limit  at  the  collector  Vx  =  0 ,  p  =  oo  ,  i  =  linjection , 

the  next  step  being  an  abrupt  drop  of  l  to  zero  as  the  magnetic  field  increases  beyond  its 
cut-off  value,  the  flow  of  electrons  across  region  II  in  Fig.  2  ceasing  altogether. 


Let  us  continue  with  our  idealised  model  of  an  infinitely  smooth  or  jelly-like  electron 
stream.  The  first  step  in  finding  point  B  is  to  note  that  the  moment  of  electron 

‘hesitancy’  (x‘  =  x"  =  0 ),  which  in  the  absence  of  the  magnetic  field  is  normally  referred 
to  as  a  ‘virtual  cathode’,  is  no  longer  associated  with  <j)  =  0 ,  but  occurs  at  some  other 
value  of  (f>  =  .  As  an  example  let  us  consider  electron  trajectories  for  6)c0  =0.5  and 

several  values  of  the  injection  current  l ,  as  shown  in  Fig.  7.  Clearly  there  is  one  value  of 
l  =  lB  for  which  an  electron  momentarily  travels  parallel  to  the  cathode,  X‘  =  0 .  We  then 
find  from  (2.9)  that  at  the  same  point  we  also  have  X“  =  0 ;  thus  just  as  was  the  case  for 
x'  =  0  and  <yc0  =  0  ,  the  electron  now  ‘hesitates’  whether  to  turn  towards  the  exit  or  the 


entrance  electrode.  Physically  this  implies  that  at  point  B  the  electric  force - EB  >0 


exactly  balances  the  magnetic  force  — <yc0y'  <0  (see  (2.3)).  In  general  point  B  can  be 
found  by  first  putting  x"  =  x'  =  0  in  (2.8),  (2.9)  and  then  expressing  Ej  and  l  in  terms 
of  rB  .  This  gives  us 


l-cB’ 


^^cOSB 

1-Cb 


(2.27) 


where 

sB  =  sin  ac0fB ,  cB  =  cos  ac0rB  (2.28) 
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(2.27)  being  a  parametric  representation  of  a  curve.  The  cross  section  of  (2.27)  and 
E,  (/)  gives  us  point  B  which  defines  the  one  particular  trajectory  for  which  momentarily 

the  x-directed  velocity  and  acceleration  disappear,  x‘  =  x"  =  0  .  For  clarity  B  points  have 
been  connected  by  a  dotted  line  in  Fig.  5.  The  corresponding  values  of  ,  <5^  ,  rd  B  ,  0B , 

E1B  and  tb  are  shown  in  Table  2.2. 

It  should  be  noted  that  for  the  cut-off  condition,  <7>c0  =  1  the  two  points  A  and  B 
coalesce,  the  large  loop  of  the  Ej  =  E,  (i)  curve  collapsing  to  a  cusp. 


2.3.3  The  new  ‘virtual  cathode’ 

In  the  presence  of  a  magnetic  field  the  conditions  for  virtual  cathode  have  to  be 
reconsidered.  We  start  by  suggesting  that  all  B  points  represent  a  highly  unstable 
solution,  which  is  similar  to  that  obtained  in  the  absence  of  a  magnetic  field.  A  more 
likely  stable  solution  is  shown  in  Fig.  8a  where,  following  the  example  of  COc0  =  0 ,  we 
have  abandoned  the  continuity  of  l  by  assuming  that  only  its  fraction  l2  reaches  the  exit 
electrode  (collector),  the  rest,  i.e.  \  =l  —  l2  turning  back  and  landing  on  the  entrance 

electrode  (anode).  We  shall  see  that  this  requires  a  shift  in  the  position  of  the  x"  =  x'  =  0 
plane  from  xB  to  some  new  position  x  ,  an  example  of  corresponding  electron 
trajectories  for  coc0  —  0.5  being  shown  in  Fig.  8b  .  However  we  will  show  that  for  a  cut¬ 
off  magnetic  field,  B  — >  BH  or  6)c0  — >  1  from  below,  the  plane  x  =  x  must  move  to  the 
exit  electrode,  making  then  x  =  x2  =  1 .  Since  the  potential  of  the  exit  electrode  is  kept 
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constant,  this  will  destroy  the  function  of  a  ‘virtual  cathode’  and  should  stop  the  vircator 
oscillations  altogether. 

When  discussing  either  a  diode  [7]  or  a  vircator  [10-12]  we  invariably  start  with  the 
assumption  that  the  beam  is  perfectly  smooth.  We  then  discover  as  we  have  pointed  out 
before  that  in  the  absence  of  the  magnetic  field  some  solutions  are  stable  and  other 
solutions  are  unstable,  the  unstable  solution  (B  in  Fig  .9)  leading  to  the  discontinuity  of 
current  flow  and  the  formation  of  a  virtual  cathode  (C  in  Fig.  9).  This  is  an  ideali2ed 
situation  based  on  an  infinitely  smooth  model  of  electron  flow.  Any  departure  from  that 
assumption  causes  instability  and  fluctuations  in  the  form  of  appearance  and 
disappearance  of  the  virtual  cathode,  thus  forming  the  basis  of  operation  of  a  vircator. 
Thus  in  practice  the  system  is  capable  of  generating  relaxation  oscillations  which  then 
generate  electromagnetic  oscillations  in  the  output  waveguide.  Fig.  1.  We  have  followed 
exactly  the  same  path  in  analyzing  the  situation  in  the  presence  of  a  magnetic  field.  In  the 
process  we  have  defined  an  equivalent  of  the  C  distribution,  but  now  in  the  presence  of 
the  magnetic  field  (Fig.  8a  &  C  in  Fig.  10b),  and  have  shown  that  the  amplitude  of  the 
corresponding  virtual  cathode  oscillations  gradually  decreases  as  the  magnetic  field 
increases. 

Let  us  now  put  the  above  brief  discussion  into  an  algebraic  form.  First  of  all  we  require 
equations  of  motion  which  would  apply  in  part  a  of  the  inter-electrode  space.  Fig.  8a. 

Since  X'  =  x"  =  0  by  definition  at  X  =  X ,  we  obtain  from  (2.7)-(2.9): 


x  = 


E, 


la 


M)+A 


s  + 


* 


cO 


co, 


cO 


2(0, 


cO 


(2.29) 
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(2.30) 


x-  =0 


_Jl 


2(0, 


■s  +  c  + 


* 


cO 


2(0, 


(1-0 


cO 


A  E,  _ _ ^  ^ 

:  0  = - -C  —  6?c0S  +  .  S 


2(0, 


cO 


(2.31) 


where  6  =  t vc0t  ,  s  =  sin«c0r  ,  c  =  cos®c0r  (2.32) 

T  being  the  electron  transit  time  from  the  entrance  electrode  at  x  =  x,  =  0  to  the  point 
x  =  x  .  Here  is  the  ‘current’  density  in  section  a  of  the  diode;  it  consists  of  the 
injection  current  &  and  the  additional  part  which  represents  the  electrons 

returning  to  the  entrance  electrode.  Therefore  from  the  point  of  view  of  electron  space 
charge  the  two  current  streams  travelling  in  opposite  directions  must  be  added  and  we 
write 


(2.33) 


since  ^  =  t2  and  <5^  =  <5^  .  In  part  b  of  the  inter-electrode  space  the  equations  of  motion 
and  the  boundary  conditions  are  somewhat  different.  Introducing  new  variables: 


X  =  x -x  ,  X'=x'-x'=x' 


Y  =  y-y,  Y‘=y -y- =oc0X 


(2.34) 
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we  can  write  in  place  of  (2,7)-(2.9): 


X  = 


2®i 


(0-S) 


(2.35) 


x=^(  1-C)  (2.36) 

2&>c0 

X"=^-S  (2.37) 

2^co 


where 

©  =  <»c0T  =  0d-0  =  ac0 (rd  -f) ,  S  =  sin©c0T ,  C  =  cosd>c0T ,  ^  ^ .  (2.38) 


We  note  that,  by  definition,  X" 


X'  =  0  in  the  plane  x  =  X  and  only  ^  — 


reaches 


the  exit  electrode.  Introducing  boundary  conditions  at  the  exit  electrode  we  now  obtain 
from  (2.35),  (2.36): 


X2  =  1  —  X 


j5_ 

2a)o0 


(®2-S2) 


(2.39) 


=  x; 


(2.40) 
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The  system  has  now  been  fully  specified,  since  there  are  five  simultaneous  equations 
(2.29)-(2.31),  (2.39),  (3.40)  and  five  unknowns:  x ,  3^  (or  ^  ),  r  ,  rd  and  E,  ,  the 


values  of  the  magnetic  field  tt>c0  and  the  injection 


<7 


At  this  point  it  is  preferable  to  reduce  the  number  of  equations  from  five  to  two.  We  do 
that  by  noting  that  by  eliminating  Ej  we  obtain  from  (2.30),  (2.31) 


^—2 

=  (2.41) 

1  —  c 

and  hence  from  either  (2.30)  or  (2.31): 


1  —  c 


Substituting  in  (2.29)  we  now  obtain 

~  <9-S 

X=~Z—? - 3T 

^co(l-c) 


(2.42) 


(2.43) 


Eliminating  x  between  (2.39)  and  (2.43)  we  have 


0 -s 
®co(l-c) 


=  l-rgr(®2-S2) 


(2.44) 
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Using  now  (2.41)  to  eliminate  S2  in  favour  of  S ,  wliich  is  given,  we  finally  obtain 


0  -s 

^co(l-c) 


®coO-c) 


(e2-s2) 


(2.45) 


Also  eliminating  S2  from  (2.39)  and  (2.40)  we  obtain,  by  substituting  for  x  from  (2.43) 


1- 


0 -s 

®co0“g) 


“2 


©0  ~S, 


CO, 


cO  _ 


CO. 


cO 


(!-c4 


(2.46) 


Using  (2.46)  the  previous  equation  (2.45)  can  be  further  simplified  by  writing 

(2-47) 


Thus  the  problem  has  now  been  reduced  to  the  solution  of  two  simultaneous  algebraic 
equations,  (2.46),  (2.47)  with  two  unknowns,  6  (or  f  )  and  02  -0A-O  (or 
„  _  4 

r2  =Td-T  ),  both  the  magnetic  field  C0c0  and  the  injection  current  S  =  —  l  being 

given.  Once  that  has  been  achieved,  the  remaining  unknowns,  i.e.  X,  (or  S2  )  and 
Ej ,  can  be  obtained  respectively  from  (2.41)-(2.43).  As  an  example  we  have  in  Table  2.3 
the  values  of  the  above  variables  for  l  =  Jmax  and  six  different  values  of  6)c0 . 
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Table  2. 1  Maximum  injection  current  at  different  values  of  the  magnetic field  &>c0 


«cO 

^max 

^nax 

Ta 

d,max 

^max 

Ei 

0.00 

8.0000 

3.5556 

1.5000 

— 

2.6667 

0.50 

7.1360 

3.1715 

1.5745 

0.78722 

2.4080 

0.80 

5.5425 

2.4633 

1.7518 

1.4014 

1.8969 

0.90 

4.6428 

2.0635 

1.8872 

1.6984 

1.5837 

0.95 

4.0082 

1.7814 

2.0071 

1.9067 

1.3491 

1.00 

2.5781 

1.1458 

2.4120 

2.4120 

0.76376 

Table  2.2  Parameters for  defining  point  B  for  different  values  of  the  magnetic field  O)c0 


^c0 

rd,B 

Em 

0.00 

4.0000 

1.7778 

3.0000 

— 

2.6667 

1.50000 

0.50 

3.5426 

1.5745 

3.1637 

0.8197 

2.3017 

1.6295 

0.80 

2.9228 

1.2990 

3.3840 

1.5562 

1.6245 

1.9344 

0.90 

2.7262 

1.2116 

3.2577 

1.9146 

1.2675 

2.1253 

0.95 

2.6651 

1.1845 

3.3065 

2.1222 

1.0620 

2.2458 

1.00 

2.5781 

1.1458 

2.4120 

2.4120 

0.7637 

2.4120 
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Table  2.3  Conditions  at  a  ‘virtual  cathode’  in  the  presence  of  a  magnetic field  a>c0 


«cO 

^max 

T 

X 

E, 

h 

0.00 

8.0000 

0.7977 

3.0000 

0.2659 

5.0145 

1.8558 

0.50 

7.1360 

0.8481 

3.2112 

0.2844 

4.6455 

1.5699 

0.80 

5.5425 

0.9724 

3.7588 

0.3308 

3.9040 

1.0718 

0.90 

4.6428 

1.0717 

4.2310 

0.3687 

3.4384 

0.8129 

0.95 

4.0082 

1.1630 

4.7025 

0.4041 

3.0822 

0.6420 

1.00 

2.5781 

2.4120 

2.4120 

1.0000 

0.7638 

(0.0000) 
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2.3.4  Numerical  examples 

Let  us  now  consider  some  numerical  examples  by  choosing  three  different  cases: 
(Dc0  =  0  (no  magnetic  field),  cdc0  =0.5  and  (Dc0  =1.0  (cut-off  condition).  In  Fig.  9  we 
have  four  well-known  potential  distributions  when  <a>c0  =0  [2]:  they  are  l  =  lmax  =8.0 
(curve  A  —  marginally  stable),  l  =  l+  where  l+  is  fractionally  greater  than  Zmax  and  thus 

generates  a  virtual  cathode  at  x=  0.2659  (see  Table  3.3)  and  the  dotted  unstable 
potential  distribution  (curve  B)  associated  with  point  B  in  Fig.  5.  This  distribution  is  by¬ 
passed  on  the  way  to  the  l+  distribution,  which  is  nominally  stable,  curve  C.  It  is  well 
known  that  if  i  =  l+  is  then  reduced,  the  virtual  cathode  persists  until  l  =  4.0  and 
x  =  0.5 ,  when  the  potential  minimum  suddenly  jumps  to  (fm  =  0.75  ,  thus  giving  rise  to  a 
new  distribution  A! .  This  generates  a  hysteresis  effect  which  was  of  great  interest  in  the 
past  [2]  and  is  now  of  fundamental  importance  in  the  operation  of  vircators. 

A  word  of  explanation  is  required  at  this  point.  According  to  our  model,  Fig.  3,  the  value 
of  the  injection  current  i  is  determined  by  the  space  charge  limited  diode  I;  for  a  fixed 
Va  we  have  1=1  when  d0  =  d  and  l  ><  1  for  d0  <>  d.  When  the  potential 

distribution  in  the  short-circuited  diode  II  follows  A  or  A'  in  Fig.  9,  then  l2  =  l  and 
hence  l  — 12  =  0 ,  no  current  flowing  from  the  anode  to  the  battery,  the  whole  of  the 
injection  current  ending  up  on  the  collector.  When  we  try  to  increase  the  injection 
current  beyond  l  =  lmax  to  what  we  have  called  l+ ,  the  system  responds  discontinuously 

and  adopts  the  potential  distribution  C.  Now  only  t2  part  of  the  ‘injection  current’  l 
reaches  the  collector,  the  rest  i.e.  l-i^  reaching  the  battery  direcdy  via  the  anode 
connection,  as  shown  in  Fig.  3.  Since  more  electrons  enter  the  short-circuited  diode  than 
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reach  the  collector,  we  allow  for  it  from  the  space-charge  point  of  view  by  writing,  using 
Figs.  8a  and  8b,  lA  =  2l  -  ^  =  l2  ■  Here  possible  misunderstanding  may  arise  since 

following  the  tradition  we  use  current  density,  which  is  a  vector  quantity,  to  simulate 
space  charge  density,  which  is  a  scalar. 

In  the  presence  of  the  magnetic  field  the  situation  is  more  complex  since  we  now  have 
(j)  =  x'2  +y'2  &  x'2  =  u2 .  Putting  <yc0  =  0.5  and  l  =  /max  =  7.1360  from  Table  2.1  we  now 
obtain  Fig.  10a,  where  again  curve  A  is  plotted  for  l  =  /max  and  curve  C  for  l  =  l+,  the 
unstable  solution  B  being  shown  dotted.  A  similar  hysteresis  effect  is  present  if  we  now 
start  reducing  i+  ,  the  zero  point  of  the  u2  -curve  creeping  along  until  it  reaches 
x  =  0.55548  ;  the  continuity  of  the  current  is  then  momentarily  restored  at 
l  =  l2  =3.5831,  but  since  the  distribution  is  unstable  the  zero  point  of  the  u2  -curve 
jumps  up  to  u2in  =  0.6800  to  form  the  minimum  of  a  stable  distribution  represented  by 
curve  A' .  In  Fig.  10b  we  have  a  similar  sequence  of  events,  but  now  expressed  in  terms 
of  the  potential  distribution  (j)  .  Here  we  note  that  due  to  the  presence  of  y  -directed 
velocities  the  potential  function  is  never  equal  to  zero,  even  when  u2  =  0 .  Also,  if  we 
compare  the  size  of  the  u2  jump  in  Fig.  10a,  it  is  only  91%  of  that  in  the  absence  of  the 
magnetic  field,  Fig.  9.  This  trend  continues  as  the  magnetic  field  increases  (e.g.  68%  for 
a>cQ  =  0.8)  and  it  is  bound  to  affect  unfavourably  the  performance  of  a  vircator. 


As  another  example  let  us  consider  the  distribution  of  the  velocity  function  u2  for 
different  values  of  the  injection  current  l  (or  .5^ )  and  a  fixed  value  of  COc0 ,  choosing 
i yc0  =  0.5  and  1.0  as  lying  at  two  ends  of  the  range  of  values  of  the  magnetic  field  for  a 


90 


conducting  diode.  We  find  that  the  value  of  the  injection  current  has  a  profound  effect 
on  the  shape  of  the  velocity  function  u2  (x) ,  as  could  be  inferred  from  Fig.  3.10.  Thus  in 
the  case  of  cdc0  =  0.5  the  curve  corresponding  to  zero  value  of  the  injection  current 
exhibits  a  negative  second  derivative,  whereas  all  the  other  curves  shown  in  Fig.  11a  have 
a  positive  second  derivative.  At  the  same  time  the  position  of  the  minimum  moves 
towards  the  centre  as  S  increases,  although  it  never  quite  reaches  it.  The  situation  is 
even  more  interesting  in  the  cut-off  case,  i.e.  for  COc0  =  1  .  Now  the  only  two  curves 
which  exhibit  a  min/max  are  those  for  1  =  0  and  l  =  the  first  having  a  maximum  at 

the  entrance  electrode,  X  =  0  and  the  second  having  a  minimum  at  the  exit  electrode, 
x  =  1 .  The  remaining  curves  exhibit  neither  a  minimum  nor  a  maximum  point,  but  only 

an  inflexion  point  at  u2  =  —  S2 ,  as  indicated  in  Fig.  3.9  for  a?c0  =  1 . 


It  is  interesting  to  note  that  unlike  the  potential  <j)  ,  some  U2  -  curves,  in  particular  those 
shown  in  Fig.  11,  exhibit  negative  curvature  over  part  of  the  inter-electrode  space.  This  is 
in  keeping  with  the  corresponding  differential  equation  which  is  obtained  from  (2.15) 
and  (2.2)  by  writing 


d2u2 

dx2 


u 


(2.48) 


The  related  Fig.  12  clearly  indicates  that  the  second  derivative  of  U2  may  have  both 
positive  and  negative  values,  the  point  of  inflexion  D  in  Fig.  1  lb  being  given  by 
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(2.49) 


_  2 1 
U  ~  2<sJ20  9ty  20 

As  the  third  example  we  must  now  consider  the  crucial  case  of  G)c0  =  1 ,  i.e.  the  ‘cut-off 
condition.  We  find  from  Table  2.3  that  for  l  =  ;max  the  point  x  has  moved  to  the  exit 
electrode,  the  corresponding  u2 -distribution  being  shown  in  Fig.  lib;  at  the  same  time 
l2  disappears  altogether.  It  is  important  from  our  point  of  view  to  note  that  x  =  x2  =  1 
and  thus  it  is  now  situated  at  the  collector  where  U2  =  0  and  the  potential  is  fixed, 
^2=1;  consequently  any  attempt  at  increasing  the  injection  current  from  zmax  to  l+ 
would  be  rejected  by  the  system,  the  excess  current  l+  —  imax  merely  flowing  back  to  the 
battery  through  the  anode  conductor.  Thus  in  this  case  the  hysteresis  effect  and  with  it 
the  generating  mechanism  for  vircator  oscillations  would  disappear  altogether. 
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2.4  Solutions  in  terms  of  distance 


In  order  to  obtain  an  alternative  description  of  our  system,  let  us  now  use  the  distance  X 
as  an  independent  variable  and  attempt  to  solve  the  corresponding  Poisson  equation.  We 

would  then  obtain  either  the  potential  <j)  or  the  velocity  x'  directly  as  functions  of  x .  In 
the  absence  of  an  initial  velocity  such  solutions  have  been  obtained  in  a  distant  past  [7,8], 
but  to  our  knowledge  they  have  never  been  extended  to  the  more  general  case  which  is 
of  interest  to  us  here,  and  which  includes  the  presence  of  a  single  initial  velocity,  in  our 

case  xj  =1. 


It  should  be  added  that  in  principle  (2.7)  —  (2.9)  are  adequate  to  obtain  all  the 
information  we  require,  since  for  any  given  r  the  distance  x  and  the  corresponding 

velocity  X'  can  be  obtained  respectively  from  (2.7)  and  (2.8),  once  0)c0  and  l  have  been 

chosen.  The  corresponding  potential  function  <j>  =<f>  (x)  can  then  be  obtained  from 

(2.15).  However,  it  seemed  desirable  for  the  sake  of  completeness  and  as  an  additional 
check  on  the  correctness  of  our  earlier  results  to  solve  the  Poisson  equation  directly.  In 
ID  the  two  basic  equation  are  given,  by  the  following  simplified  expressions 


d V 

dx2 


dE  _  . 

- ~-p>  0 

dx 


Poisson  equation 


(2.50) 


—  Jc  =  —  (/?X‘)  =  0  Continuity  equation  (2-51) 

dx  dx 


The  latter  equation  establishes  the  following  important  property  of  Jc  in  ID 
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Jc  =  const 


(2.52) 


Substituting  (2.15)  in  the  Poisson  equation  (2.50)  and  using  our  definition  for  the  current 
density  (2.5)  we  obtain 


ay  _ 

dx2  x  u 


s 


(2.53) 


Here  we  have  used  S'  —  —  Jc  and  u  =  x'  for  typographical  convenience.  At  this  point  it 

is  important  to  note  that  in  the  absence  of  the  magnetic  field  U2  =  <j)  ,  so  that  the  two 
functions  are  identical.  This  is  no  longer  the  case  in  the  presence  of  the  magnetic  field,  as 
is  indicated  by  (2.15).  Since  now  the  appearance  of  a  virtual  cathode  is  associated  with 
u2  =0  and  not  <j)  =0 ,  we  must  alter  (2.53)  in  order  to  obtain  u2  =  u2(x)  rather  then 

</)=<!>  (x)  .  This  also  conveniently  removes  the  independent  variable  x  from  the 
denominator  on  the  right  hand  side  of  (2.53).  Differentiating  (2.15)  twice  with  respect  to 
X  we  obtain. 


dV 

dx2 


d2u2 

dx2 


+2ty20 


(2.54) 


and  after  substitution  in  (2.53) 


d2u2 

dx2 


u 


(2.55) 
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At  this  point  it  is  important  to  note  that  the  transition  from  (3.54)  to  (3.55)  i.e.  from  (/) 
to  u2  is  by  no  means  trivial.  Comparing  Fig.  12  and  Fig.  4  we  find  that  although  the 
second  derivative  of  (f>  must  be  always  positive  in  the  presence  of  negatively  charged 
electrons,  this  no  longer  applies  to  the  second  derivative  of  u2  when  the  magnetic  field  is 


du2 


present,  although  there  is  a  limit  on  its  negative  values  given  by  — » -2(2)20  as 

dx 


U  —>  oo  ,  the  point  of  inflection  occurring  when  u  =  — —  .  All  this  is  direcdy  related  to 

2®c0 

the  fact  that  in  the  presence  of  a  magnetic  field  we  have  at  any  given  point  two 

1  -  .  _ _2  _ 

competing  forces,  the  electric  force  — — E  and  the  magnetic  force  a>c0y' =  O)c0x  ,  as 
indicated  in  (2.3). 


Multiplying  both  sides  of  (2.55)  by  du2/dx  and  integrating  with  respect  to  X  we  obtain 


1  du2  _du  rzn  zyzi 

=  u— =  V^  u-t»c0u  +C,  (2.56) 

2  dx  dx 

It  is  interesting  to  compare  (2.56)  and  (A2.15);  they  are  identical  except  for  the  presence 
of  an  extra  term  — ty20 U2  due  to  the  presence  of  the  magnetic  field.  It  should  be  noted 

that  this  term  opposes  the  effect  of  space  charge  expressed  by  the  S  U  term.  A  similar 
term  appears  in  the  presence  of  positive  ions  [13],  so  that,  as  has  already  been  pointed 
out  in  Section  2.2,  the  presence  of  the  magnetic  field  has  the  effect  of  increasing  the 
‘perveance’  of  the  system. 

A  second  integration  with  respect  to  X  now  leads  to  the  following  expression 
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(2.57) 


“W^  u-^cou'  +  Ci-^-s111 

^cO  2®c0 


2ft>c0u-^ 

V^+^ic, 


=  x  +  C2 


The  above  equation  substantially  differs  from  (A2.16),  largely  due  to  the  presence  of  the 


—2 


sin  1  function,  which  makes  it  possible  for  — to  become  negative  for  u  >  y / 2co20  , 

dx 

as  shown  in  Fig.  3.10.  This  is  quite  a  marked  departure  from  the  case  of  B  =  0  and  is  the 
basis  of  the  so-called  striations  observed  by  others  [4]. 


The  first  constant  of  integration  Q  can  be  easily  eliminated  by  nothing  that  U  =  1  when 
X  =  0 .  Hence  (2.57)  reduces  to: 


=4rU^u-(0 cou'+C,  -yl<?-G)>0+cX 


y 

2®c0 


2o)2c0u-y 

V^2+4«c20C1 


■sin 


2<^c0  ^  > 

\]y2  +  46>20c, 


(2.58) 


The  other  constant  of  integration  is  more  difficult  to  eliminate.  In  fact,  following  an 
earlier  example  [2]  indicated  in  (A2.19),  we  can  choose  a  new  constant  which  is  Um ,  i.e. 
the  min/max  value  of  the  velocity  u,  when  it  exists  between  the  electrodes.  Writing 
du2/dx  =  0  on  the  LHS  of  (2.56)  we  obtain 


o  =  y  u„ 


■<y2”2 


c0 


u»+Ci 


(2.59) 


or 
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(2.60) 


u„  =■ 


2fiic0L 


2+4<o 


ci 


Substituting  in  (2.57)  we  obtain 


=  4-U^(u-um)-tyc20(u2-<)-^(l-um)-(yc20(l-u^)j 

‘"cO  ^  > 


3T 


sin 


-i  2®, 


"2\i~y  . 


c0‘ 


2®cOUm-^ 


-sin 


2«c2o-^ 


(2.61) 


2«c20um 


One  might  argue  that  in  the  presence  of  a  min/max  between  the  electrodes  the  effect  of 
the  constant  of  integration  is  somewhat  easier  to  comprehend  in  (2.61)  than  (2.58). 

Let  us  now  consider  some  examples  of  the  potential  and  velocity  distribution  between 
the  electrodes.  In  the  absence  of  a  magnetic  field  we  have  the  well  known  distributions 
shown  in  Fig.  A4.  Since  in  this  case  (/)  =  u2 ,  the  two  sets  of  distributions  are  identical; 
they  are  all  symmetrical  with  a  minimum  situated  exactly  in  the  middle  plane  at  X  =  0.5 . 
The  situation  is  somewhat  different  when  the  magnetic  field  is  present.  In  Figs.  3.11  and 
3.12  we  have  respectively  the  two  function  (f)  (x)  and  U2(x)  for  l  -  tmax  and  different 
values  of  <2>c0  .  We  find  from  Fig.  6a  that  as  d)c0  increases  from  zero  to  1.0  (cut-off 
condition),  the  position  of  the  potential  minimum  xm  shifts  from  the  center  towards  the 
exit  electrode,  however  the  potential  at  the  entrance  electrode  reaches  unity  for  all  values 
of  <yc0 .  The  situation  is  dramatically  different  in  the  case  of  velocity  distribution.  Due  to 

the  presence  of  the  magnetic  field  the  value  of  u2  at  the  exit  electrode  is  given  by 
1  -  20 ,  as  shown  in  (2.16);  in  fact  it  is  zero  for  a>c0  =1 ,  i.e.  for  the  cut-off  condition. 
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Also  the  shift  of  the  minimum  of  the  U2  function  away  from  the  middle  plane  is  much 
greater  than  in  the  case  of  the  potential  function.  In  fact  it  can  be  shown  that  for  (Dc0  =  1 

the  position  of  the  potential  minimum  of  the  U2  function  coincides  with  the  exit 
electrode,  so  that  Xm  =  X2  =  1 . 

As  another  example  let  us  consider  the  distribution  of  the  velocity  function  u2  for 
different  values  of  the  injection  current  i  (or  S )  and  a  fixed  value  of  <0cO ,  choosing 

a>c0  =  0.5  and  1.0  as  lying  at  two  ends  of  the  range  of  values  of  the  magnetic  field  for  a 
conducting  diode.  We  find  that  the  value  of  the  injection  current  has  a  profound  effect 
on  the  shape  of  the  velocity  function  u2  (x) ,  as  could  be  inferred  from  Fig.  12.  Thus  in 
the  case  of  6)c0  =0.5  the  curve  corresponding  to  zero  value  of  the  injection  current 
exhibits  a  negative  second  derivative,  whereas  all  the  other  curves  shown  in  Fig.  6b  have 
a  positive  second  derivative.  At  the  same  time  the  position  of  the  minimum  moves 
towards  the  centre  as  S  increases,  although  it  never  quite  reaches  it.  The  situation  is 
even  more  interesting  in  the  cut-off  case,  i.e.  for  G)c0  =  1  .  Now  the  only  two  curves 
which  exhibit  a  min/max  are  those  for  i  =  0  and  l  =  tmax  ,  the  first  having  a  maximum  at 
the  entrance  electrode,  X  =  0  and  the  second  having  a  minimum  at  the  exit  electrode, 
X  =  1 .  The  remaining  curves  exhibit  neither  a  minimum  nor  a  maximum  point,  but  only 

an  inflexion  point  at  u2  =  —  S2 ,  as  indicated  in  Fig.  8a  for  <ac0  =  1 . 

It  should  be  added  that  most  of  the  curves  shown  in  Figs.  6a,  b  and  11  have  been 
obtained  using  time  as  the  independent  variable,  as  explained  in  Section  2.3,  although 
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they  then  have  been  checked  against  the  much  more  cumbersome  expressions  derived  in 
the  present  Section. 


2.4.1  Zero  magnetic  field 

Let  us  start  our  investigation  with  the  simpler  case  of  a  zero  magnetic  field,  B-  0 .  The 
question  of  stability  of  electron  flow  in  a  diode,  when  the  system  is  not  space-charge- 
limited,  came  into  prominence  in  the  late  thirties  in  connection  with  the  development  of 
multi-grid  tubes,  the  power  tetrode  in  particular.  The  problem  was  formulated  in  two 
famous  papers  by  Llewellin  and  Llewellin  &  Peterson  [18],  where  the  steady  state 
solutions  given  by  Fay  et  al.  [2]  have  been  perturbed  by  a  small  a.c.  injection  current.  The 
relationship  between  the  injection  conditions,  and  the  corresponding  output  potential, 
current  and  velocity  (all  ac)  is  expressed  by  the  well-known  3x3  matrix  of  ai}  coefficients. 

An  excellent  summary  of  the  derivation  of  the  atj  coefficients  is  provided  in  Section  3b 
of  [19].  A  more  detailed  discussion  of  the  stability  of  such  solutions  based  on  the  energy 
considerations  and  also  on  the  algebraic  properties  of  one  of  the  a y  coefficients,  viz.  an 

is  provided  in  Ch.  3  of  [7].  As  mentioned  before  a  physically  convincing  ‘explanation’  of 
the  source  of  possible  instability  of  the  system,  using  a  minimum  of  mathematics,  in 
Section  10.6  of  [8]. 

The  an  coefficient  of  the  atJ  matrix  relates  the  dc  injection  current  J c  to  the  ac 
potential  difference  between  electrodes  <j>2~^> i  •  Using  our  bar  notation  we  can  write: 


tj)2  <f>x  -  auJc 


(2.63) 
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where  from  [19,  p.  37] 


fli  i  — 


Zd  L \+V2 


10) 


—  __d_ 

id) 


l-lf(l  +  2  ®() 


1 


l--C(l  +  20>4) 


(2.64) 


since  for  a  short-circuited  diode  we  have  in  our  bar  notation  (see  Appendix  1)  the 
injection  velocity  Dj  =  u2  =  1  .  Here  fd  is  the  steady  state  transit  time  between  the 
electrodes  and  C,  is  the  so-called  space-charge  factor  given  by  two  alternative  definitions: 


(  f  > 

1  VO 

=  3 

U-4-) 

<  *d  J 

l  TJ 

(2.65a) 


„  1  - 

C  = - —  =  —  IT 

4(t71+t72)  18 


(2.65b) 


where  Td0  is  the  value  of  Td  in  the  absence  of  space  charge  and  l  =  ——Jc  is  the 


normalized  injection  current  density,  as  defined  in  Appendix  1  and  then  used  in  Fig.  5. 
Substituting  (2.65b)  in  (2.64)  we  obtain 


°ii 


ICO 


_ 

id) 


1  2 

'-3^3^ 
fj  3  18  * 


—3 


IT. 


llico 


27  ^ 


(2.66) 
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Since  for  a  shortcircuited  diode  we  have  <j)2  =  ^  ,  the  1-h.s.  of  (2.63)  is  zero  and  we  thus 
obtain  from  (2.66)  the  required  relationship 


=0  or  A,  =  t 1>4, 


(2.67) 


where  <J>4  is  given  by 

04=-4(2-z-2e"z-ze"2)  (2.68) 

Thus  the  question  of  stability  of  the  solutions  has  now  been  reduced  to  the  investigation 
of  (2.67).  In  the  initial  derivation  of  (2.64)  the  variable  z  was  assumed  to  be  purely 
imaginary,  z  =  iCDt ,  since  it  related  direcdy  to  the  initial  sinusoidal  perturbation  of  the 
system  by  adding  an  ac  term  to  the  dc  Llewellyn  -  Peterson  equation  (2.6).  Taking  z  to  be 
pure  real  and  introducing  some  approximations  to  <X>4  the  stability  of  solutions  has 
already  been  discussed  by  Birdsall  [7].  We  decided  to  avoid  any  approximations  to  <1>4  as 
being  too  dangerous,  bearing  in  mind  that  in  the  presence  of  a  magnetic  field  the 
equivalent  of  ti>4  acquires  quite  a  complicated  form,  (2.76).  We  have  therefore  adopted  a 
more  general  approach  by  treating  z  as  a  complex  variable  z  =  x  +  iy ,  keeping  (2.68)  as  it 
stands  and  simply  looking  for  zeros  of  (2.67). 

The  real  and  imaginary  parts  of  <I>4  are  respectively  shown  in  Figs.  13  a  and  b.  We  find 
that  in  spite  of  its  apparent  algebraic  complexity  the  function  can  be  envisaged  as  a 
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simple  combination  of  exponentials  along  the  real  x-axis  and  trigonometric  oscillations  in 
the  imaginary  iy  -direction. 

The  solution  of  (2.67)  is  simplified  by  the  fact  that  the  l.h.s.  is  pure  real,  the  values  of  l 
and  rj  directly  relating  to  Fig.  5.  Thus  in  order  to  solve  (2.67)  we  simply  plot  the 
contours  Re  <I>4  =  and  Im  <F4  =  0 ,  the  points  where  the  two  sets  of  contours  cross 
being  the  required  solution. 

As  an  example  let  us  consider  point  A  in  Fig.  5,  which  represents  the  maximum 
permissible  value  of  the  injection  current  in  a  short-circuited  diode  given  by 
1  =  lmm  =  8.0  ,  Table  2.1.  For  CDco  =  0 ,  (2.21)  acquires  a  much  simpler  form  given  by  (see 
A2.8). 

(2-®) 

rd 

Solving  (2.69)  for  l  =8.0  we  obtain  Td  =1.5 ;  substituting  for  both  l  and  Td  in  (2.55) 
we  find  that 

A,=~  =  1.0000  =  0>4  (2.70) 

The  superposition  of  the  two  sets  of  contours  Re  <F4  =1.0,  Im  <I>4  =0.0  is  shown  in 
Fig.  14a,  the  corresponding  positions  of  the  zeros  being  indicated  by  circles.  As  was  to  be 
expected  the  important  main  zero  lies  at  the  origin  of  the  complex  z-plane,  indicating 
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that  point  A  represents  a  conditional  equilibrium  separating  the  stable  and  unstable 
regimes  of  the  system. 

In  order  to  familiarize  ourselves  with  our  technique  of  testing  the  stability  of  solutions, 
let  us  briefly  consider  two  more  points,  one  stable  (point  A’)  and  one  unstable  (point  B). 
We  find  from  (A2.8),  Table  2.2  and  Fig.  5  that  the  important  parameters  for  the  two 
points  are  given  by: 


Table  2.4  Parameters  of  Points  A’  and  B 


i 

E, 

4 

A’ 

1.0981 

4.0000 

0.9762 

5.0977 

B 

3.0000 

4.0000 

2.6667 

0.2500 

We  can  now  repeat  our  procedure  for  solving  (2.67)  with  new  values  of  .  Figs.  14b 
and  c  respectively  show  the  contours  of  Re  <X>4  for  the  two  values  of  given  in  Table 

2.4  and  for  Im  ®4  =  0.0 .  Again  the  roots  of  (2.67)  occur  where  the  two  sets  of  contours 
cross.  We  now  find,  as  to  be  expected,  that  for  point  A’  the  main  root  of  (2.67)  lies  on 
the  negative  part  of  the  real  axis,  indicating  that  the  solution  is  stable.  Similarly,  the  main 
root  of  (2.67)  for  point  B  lies  on  the  positive  part  of  real  axis,  indicating  that  the  solution 
is  unstable. 

At  this  point  it  should  be  clarified  that  the  above  analysis,  although  useful,  is  highly 
idealized.  In  practice  an  electron  stream  is  not  a  homogeneous  jelly,  but  consists  of  a 
large  number  of  individual  particles  which  contribute  to  noise,  the  other  contributions 
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being  random  velocities  and  random  times  of  emission  of  individual  electrons  [20].  This 
means  that  in  practice  the  unstable  branch  A-B  can  never  be  reached,  the  system 
breaking  out  into  relaxation  oscillations,  forming  the  basis  of  operation  of  a  vircator  [7]. 
However  the  idealized  solutions  based  on  (2.6)  and  Llewellyn  -  Peterson  equations  are 
useful  in  proving  the  fundamental  instability  of  physical  situations  they  represent.  Similar 
comments  apply  when  the  magnetic  field  is  present. 


2.4.2  Magnetic  field  present 

In  the  previous  section  we  have  considered  the  case  of  zero  magnetic  field  as  a  way  of 
introduction  to  our  relatively  straightforward  method  of  investigating  the  stability  of 
solutions.  Let  us  now  consider  the  more  difficult  case  when  the  magnetic  field  is  present. 
Fortunately  considerable  amount  of  work  has  been  done  on  the  extension  of  the 
Llewellyn-Peterson  equations  to  the  case  of  crossed-field  geometries  and  we  are  going  to 
rely  heavily  on  the  results  of  those  investigations  [21,  22].  The  main  difference  is  that  the 
magnetic  field  forces  us  to  consider  the  current  flow  in  two  dimensions.  Consequentiy 
the  usual  perturbation  of  the  steady  state  equation  (2.6)  now  leads  to  a  set  of  atj 

coefficients  which  form  a  4x4  rather  than  3x3  matrix,  as  they  did  before.  Since  there  are 
no  a.c.  components  of  current  or  velocity  at  the  entrance  electrode,  the  d.c.  injection 
current  being  perpendicular  to  the  entrance  electrode,  (2.61)  still  holds,  except  that  now 
is  much  more  elaborate  than  the  corresponding  expression  in  (2.64).  Using  our  bar 
notation  we  obtain  from  [21,  p.  82]: 
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(2.71) 


where 


P  =  \-i~sgce-i6-cece-ie 
co„ 


Q  =  l-i^sdce-i0-cece-m 

CO 


(2.72) 


and 


6  =  cot d,dc=  coc0rd  ,  s0c  =  sin ac0Td  ,  c0c  =  cos coc0Td  (2.73) 

The  new  expression  for  an  looks  quite  formidable,  fortunately  it  can  be  reduced  to  a 
more  manageable  form  for  our  particular  task  of  stability  investigations. 


In  (2.71)  the  coefficient  an  is  expressed  in  terms  of  the  initial  velocity  x{  and  initial 

acceleration  xj"  .  It  is  more  convenient  for  our  purpose  to  express  it  in  terms  of 
interelectrode  distance  xd  which  in  the  bar  notation  is  equal  to  unity  by  definition  (see 
(Al.l)  in  Appendix  1).  Thus  we  can  combine  the  first  four  terms  on  the  r.h.s.  of  (2.71) 
and  write  with  the  help  of  (2.7): 
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where  Xd  =  1  .  Now  the  last  two  terms  of  (2.74)  can  be  combined  together.  Putting 


— -  4  2 

Jc  =  -—l,  as  indicated  in  (A1.2)  and  writing  -d)2Td  =  {io)Td)  =  z2  we  obtain  from 
(2.74) 
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where 
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(2.76) 


The  main  difference  between  (2.76)  and  (2.68)  is  that  <P4e  contains  6C  =  O)c0Td  ,  where 
d)c0  refers  to  the  strength  the  magnetic  field  as  defined  in  (Al.l)  of  Appendix  1.  Figs. 
17a  and  b  respectively  show  the  real  and  imaginary  part  of  <t>4c  for  O)c0  =0.8 .  Strangely 
enough  in  spite  of  the  apparent  algebraic  complexity  of  (2.76),  the  difference  in  the 
appearance  of  the  corresponding  surfaces  in  Figs.  17a,  b  and  13a,  b  is  not  that  great.  It 
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should  be  added  that,  as  was  to  be  expected  04c  — >  ®4  as  0C  -»  0  .  Since  the  crossed- 
field  diode  we  are  considering  is  short  circuited,  as  indicated  in  Fig.  3,  the  l.h.s.  of  (2.63) 
is  again  zero,  leading  to 

3,,=o  ot4=a>4t.4=-^-  (2.77) 

i*d 

which  is  the  equivalent  of  (2.67)  when  the  magnetic  field  is  present. 

The  real  and  imaginary  parts  of  <D4c  are  respectively  shown  in  Figs.  15a  and  b  for 
(Dcq  =  0.8 .  It  is  interesting  to  note  that  in  spite  of  the  algebraic  diversity  of  (2.68)  and 
(2.76),  corresponding  surfaces  representing  their  real  and  imaginary  parts  are  remarkably 
similar,  as  can  be  seen  from  Figs.  13a,  b  and  15a,  b. 

In  order  to  obtain  the  roots  of  (2.77)  for  the  three  points  A,  A’  and  B  and  for  different 
values  of  (Dc0  ,  we  follow  the  same  procedure  as  before,  except  that  instead  of  using  <D4 
we  now  use  as  indicated  in  (2.77),  the  corresponding  values  of  l  and  Td  being 
obtained  from  (2.21),  (2.22)  and  Tables  2.1  and  2.2.  The  results  of  our  computations  are 
shown  respectively  in  Figs.  16a-c.  In  Fig.  16a  we  have  the  roots  of  (2.76)  for  point  A  and 
for  six  different  values  of  the  magnetic  field,  including  d)c0  =  0  .  We  find  that  the  main 
root  firmly  remains  at  the  origin,  the  solution  for  i  =  zmax  being  conditionally  stable, 
irrespective  of  the  corresponding  values  of  the  magnetic  field. 

In  the  case  of  point  A’  which  is  situated  on  the  lower  branch  of  the  curves  shown  in  Fig. 
5,  the  main  root  is  always  placed  on  the  negative  part  of  the  x-axis,  as  shown  in  Fig.  16b, 
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whatever  the  value  of  the  magnetic  field.  The  only  exception  is  the  limiting  case  6)c0  =  1 , 
when  the  three  points  coalesce,  the  corresponding  root  appearing  at  the  origin. 

Finally  the  main  roots  of  (2.77)  for  point  B  are  all  placed  on  the  positive  part  of  real  axis 
whatever  the  value  of  the  magnetic  field,  the  corresponding  solutions  being 
fundamentally  unstable.  The  only  exception  again  is  the  ‘cut-off  condition  when  (5c0  =  1 ; 
now  the  A’-A-B  loops  of  Fig.  5  collapse  to  a  cusp,  the  solution  becoming  singular.  It 
should  be  added  that  the  A-B  branch  of  each  curve  is  the  equivalent  of  the  well-known 
‘overlap’  solution  fully  discussed  for  a)c0  =  0  elsewhere  [2,7]. 

The  above  results  extend  the  stability'  investigations  to  a  conducting  crossed-field  diode 
and  thus  confirm  the  validity  of  the  assumptions  made  in  a  recent  publication  [23,  special 
issue  on  High-Power  Microwave  Generation]  which  covered  the  results  reported  in 
Section  2.3.  The  stability  investigations.  Section  2.5,  have  been  presented  in  a  poster 
session  at  the  IEEE  Int.  Conf.  IVEC2002  [24].  It  should  be  added  that  the  method  of 
stability  investigations  described  above  is  more  general  than  it  would  appear  from  (2.67) 
or  (2.77),  since  it  is  not  limited  to  a  short-circuited  diode.  In  the  presence  of  a  load 

au  ^  0  and  (2.63)  can  then  be  suitably  adjusted  to  allow  for  <j2  —  <f)x  ^  0 . 
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2.5  Stability  considerations 

In  Section  2.3  we  have  assumed  that  some  solution  of  (2.6)  are  stable  or  conditionally 
stable  and  that  other  solutions  are  not.  We  will  now  attempt  to  prove  that  this  indeed  is 
the  case.  It  is  necessary  for  us  to  do  so,  because  the  only  existing  proofs  apply  to  the 
simpler  case  when  the  magnetic  field  is  not  present  [7,  8]. 

The  problem  of  stability,  or  lack  of  it,  has  a  long  history.  It  all  started  with  the  question 
of  stability  of  planetary  orbits,  possibly  first  raised  by  Newton,  and  in  greater  detail  by 
Kepler.  It  was  then  defined  with  great  precision  by  Poineare  in  his  magnum  opus  on 
celestial  mechanics  [14].  Poineare  also  indicated  how  to  treat  stability  in  more  general 
systems  than  those  based  on  particle  orbits.  Almost  at  the  same  time  Liapunov  wrote  his 
seminal  thesis,  now  available  in  English  [15],  on  the  general  approach  to  the  question  of 
stability  in  a  set  of  equations  of  the  form 

=  . n  (2.62) 

It  is  interesting  to  note  that  in  the  work  on  chaos,  where  the  question  of  stability  is  of 
crucial  importance,  it  is  Liapunov  approach  in  the  form  of  Liapunov  exponentials  which 
seems  to  prevail  [16].  It  may  not  be  out  of  place  to  repeat  here  a  quotation  from  one  of 
the  best  textbooks  on  classical  mechanics  [17]:  ‘the  general  subject  of  stability  of  motion 
[of  a  system]  is  today  a  very  active  field  of  investigation,  often  using  sophisticated 
mathematical  tools....’  I  am  sure  the  author  must  have  had  TAKOMAK  and  plasma 
physics  in  mind  when  he  was  writing  it.  By  contrast,  in  our  investigations  we  are  going  to 
use  the  simplest  mathematics  available.  In  many  cases  the  problem  of  stability  can  be 
formulated  as  a  small  perturbation  of  the  system  which  may  well  take  the  form  of  a 
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dispersion  equation.  The  analysis  of  the  (complex)  roots  of  the  perturbed  system  then 
gives  us  some  idea  whether  the  system  is  stable  or  unstable.  It  is  this  approach  we  are 
going  to  adopt  in  our  investigations. 
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Conclusions  -  part  2 


A  plane  magnetic  diode  has  been  used  as  a  simplified  model  of  the  effect  that  an  axial 
magnetic  field  would  have  on  the  operation  of  an  axial  vircator.  Having  discussed  in 
some  detail  the  electron  trajectories,  we  were  able  to  indicate  a  hitherto  overlooked 
solution  reminiscent  of  a  virtual-cathode  distribution  common  in  the  absence  of  a 
magnetic  field.  However,  due  to  the  presence  of  the  magnetic  field,  u2  =  0  does  not  imply 
(ft  =  0  and  the  ‘vircator’  effect  is  weakened.  In  fact  when  the  magnetic  field  reaches  its 
near  cut-off  value,  ooc0  -»  0 ,  the  point  u2  =  0  moves  to  the  exit  electrode,  where  the 
potential  is  fixed,  <j)  =  1 .  This  must  suppress  vircator  oscillations  altogether,  an  effect 
which  has  been  observed  in  practice.  A  further  increase  of  the  magnetic  field  beyond  its 
cut-off  value  would  stop  the  current  flow  altogether. 
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Figures  -  part  2 


Figure  1  Cross-section  of  a  typical  coaxial  vircator. 

Figure  2  Radial  electron  velocity  vr  in  the  middle  plane  of  the  anode  as  a  function  of 
radius  r. 

Figure  3  The  model  consisting  of  a  space-charge  limited  diode  I  followed  by  a  short- 
circuited  crossed-field  diode  II. 

Figure  4  Poisson  equation  without  (1)  and  with  (2)  the  magnetic  field  B  expressed  in 
terms  of  (Dc0 . 

Figure  5  The  electric  field  at  the  entrance  electrode  Ej  as  a  function  of  the  injection 
current  l  for  different  values  of  the  magnetic  field  expressed  by  (Oc0  . 

Figure  6a  Potential  distribution  </)  =  <j)  (x)  for  l  =  lmsx  and  different  values  of  (Oc0  . 

Figure  6b  Velocity  distribution  u2  =  u2  (x)  for  l  =  lmsx  and  different  values  of  O)c0  . 

Figure  7  Electron  trajectories  for  (Oc0  =  0.5  and  several  values  of  the  injection  current  l . 

Figure  8a  A  possible  distribution  of  current  for  i  >  tmax  . 

Figure  8b  Electron  trajectories  for  d)c0  =  0.5  and  three  different  values  of  the  injection 
current  l . 

Figure  9  Four  critical  (f>  =  u2  distributions  for  6)c0  =  0  . 

Figure  10a  Four  critical  velocity  distributions  for  a ?c0  =  0.5 . 

Figure  10b  Four  critical  potential  distributions  for  d)c0  =0.5 . 

Figure  11a  The  dependence  of  the  velocity  function  u2(x)  on  the  injection  current  i 
for  6)c0  =  0.5 . 

Figure  lib  The  dependence  of  the  velocity  function  U2(x)  on  the  injection  current  l 
for  =1.0. 

Figure  12  Following  (2.48):  the  second  derivative  d2u2/dx2  as  a  function  of  U  . 

Figure  13a  Real  part  of  04  as  a  function  of  (x,y). 

Figure  13b  Imaginary  part  of  <J>4  as  a  function  of  (x,y). 
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Figure  14a  Contours  of  Re  <X>4  =  Ay  (full  line)  and  Im  04  =  0  (broken  line)  for  point  A; 
roots  of  (2.67)  are  marked  by  circles. 

Figure  14b  Contours  of  Re  <I>4  =  Ay  (full  line)  and  Im  <J>4  =  0  (broken  line)  for  point  A’; 
roots  of  (2.67)  are  marked  by  circles. 

Figure  14c  Contours  of  Re  <t>4  =  Ay  (full  line)  and  Im  ®4=0  (broken  line)  for  point  B; 
roots  of  (2.67)  are  marked  by  circles. 

Figure  15a  Real  part  of  04c  as  a  function  of  (x,y)  for  (0cQ  =  0.8  . 

Figure  15b  Imaginary  part  of  C>4c  as  a  function  of  (x,y)  for  G>c0  =  0.8 . 

Figure  16a  Point  A:  roots  of  (2.76)  for  six  different  values  of  the  magnetic  field 
(conditionally  stable  solutions). 

Figure  16b  Point  A’:  roots  of  (2.76)  for  six  different  values  of  the  magnetic  field  (stable 
solutions). 

Figure  16c  Point  B:  roots  of  (2.76)  for  six  different  values  of  the  magnetic  field  (unstable 
solutions). 
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Figure  3  The  model  consisting  of  a  space-charge  limited  diode  I  followed  by  a  short-circuited  crossed-field 
diode  II. 
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Figure  5  The  electric field  at  the  entrance  electrode  E!  as  a  function  of  the  infection  current  i  for  different 


values  of  the  magnetic field  expressed  bj  (Dc0  . 
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Figure  13a  Real part  of  <t>4  as  a  function  of  (x,y). 
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Figure  13b  Imaginary  pari  of  <t>4  as  a  function  of  (x,y). 
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Figure  14b  Contours  of  Re  04  =  Aj  (full  line)  and  Im  <f>4  =  0  (broken  line)  for  point  A';  roots  of 
(2.67)  are  marked  by  circles. 
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Figure  14c  Contours  of  Re  <1 >4  =A(I  (full  line)  and  Im  <D4  =  0  (broken  line)  for  point  B;  roots  of 
(2.67)  are  marked  by  circles. 
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Figure  15a  Real part  of  <X>4c  as  a function  of  (xj)  for  a>eQ  =  0.8 . 
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Figure  15b  Imaginary  part  of  04c  as  a  function  of  (x,y)  for  O)c0  =  0.8 . 
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Future  work 


One  area  which  in  our  opinion  requires  further  investigation  is  the  mode  structure  of  the 
microwave  output.  Since  the  output  waveguide  is  heavily  overmoded  and  the  frequency 
spectrum  of  the  microwave  pulse  frequently  shows  the  presence  of  additional 
resonances,  a  series  of  careful  measurements  using  suitable  waveguide  probes  would  be 
required. 

As  far  as  the  more  general  area  of  improving  the  performance  of  the  vircator  is 
concerned,  then  we  have  two  suggestions  to  make: 

1.  In  our  opinion  it  is  essential  to  pursue  the  investigations  of  the  electron  cloud 
formation  in  the  important  interaction  region  where  the  virtual  cathode  is  being 
formed.  This  task  is  substantially  facilitated  by  the  technique  of  ‘moving  pictures’ 
as  developed  by  one  of  us  (X.  Chen)  and  described  in  the  report. 

2.  Associated  with  this  we  feel  it  is  necessary  to  investigate  even  more  closely  the 
design  of  the  space  still  inside  the  anode  drum,  but  nearer  to  the  flange.  It  is  there 
that  a  cleverly  designed  septum  should  catch  stray  electrons  and  thus  improve  the 
power  output  and  electronic  efficiency  of  the  device.  All  this  would  invariably 
have  a  beneficial  effect  on  improving  the  appearance  of  the  frequency  spectrum. 
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Appendix  1.  ‘Bar’  notation 


In  order  to  reduce  the  number  of  independent  variables  it  is  convenient  to  use  the 
following  notation  based  on  some  earlier  work  [1]. 

x  =  x/d,  T  =  t/tp  v  =  v/Vj,  <f>  =  <f>l<f>x, 
p  =  (d2  le0<f>x)p,  E  =  (d/^)E  =  -5x/5^ 

Jc=pv  =  (d2 /£jxvx)Jc  < 0,  ac0  =  txo)c0 

Here 

d 

tj  =  d/v, 

®  co  =(e/m0)B 

B 

Thus  following  (A  1.1)  we  always  have: 

<f>{=  1,  V]  =(dx/dt)j  =xj  =1  at  the  entrance  electrode,  Xj  =0  and  x2  =  1  at  the  exit 
electrode.  In  addition  for  a  short-circuited  diode  <f>2  =  1  at  the  exit  electrode.  It  is  also 
convenient  to  express  the  injection  current  in  terms  of  a  reduced  current  ‘iota5,  where 


separation  of  the  electrodes 
electron  transit  time  in  the  absence  of  fields 
electron  cyclotron  (angular)  frequency 
axial  magnetic  field 


(Al.l) 


and 


(A1.3) 


is  the  current  density  in  an  equivalent  space-charged  limited  diode  [1]. 


For  purely  algebraic  reasons  we  have  also  introduced  a  ‘curly’  capital  I  defined  as 


(A1.4) 


This  simplifies  the  notation  in  some  cases. 
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Appendix  2.  Summary  of  results  for  a  short  circuited  diode  when  B=0 


As  a  check  on  the  analytical  expressions  derived  in  the  presence  of  a  magnetic  field,  B^O, 
it  is  desirable  to  compare  them  with  similar,  but  well-known  expressions,  for  the  case  of 
B^O  [1].  The  corresponding  expressions  are  summarised  here  for  the  convenience  of  the 
reader.  As  is  well  known  in  the  absence  of  a  magnetic  field  the  equation  of  motion  and 
the  energy  balance  equation  acquire  a  simplified  form;  thus  putting  a>c0  =  0  in  (3.2.1)  and 
(3.2.13)  we  obtain 


Equation  of  motion 


(A2.1) 


2  _ 

X-  =</>  Energy  balance  (A2.2) 

However  the  Poisson  equation  (3.3.1)  and  the  continuity  equation  (3.3.2)  remain 
unaltered. 


A2.1  Time-dependent  solutions 


Let  us  repeat  the  definition  of  the  current  density  (3.2.3)  by  writing 


3c  =  pic=  —  x-=E-  (A2.3) 

dx 


Differentiating  (A2.1)  with  respect  to  t  and  substituting  from  (A2.3)  we  obtain 
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x‘"  =-— E‘  =  J.  =  const 

2  2  c 


(A2.4) 


which  is  the  steady-state  part  of  the  Llewellyn-Peterson  equation  [2].  The  solution  of 
(A2.4)  is  quite  simple  and  is  given  by 


x  -t  E,r2  +— r3 
4  1  27 


(A2.5) 


x’  =1-— E,r  +—t2 

2  9 


(A2.6) 


-  ^2l~ 

x  = — E,  + — T 

2  9 


(A2.7) 


_  (, 

Here  we  used  the  same  variables  as  in  (3.2.5)  —  (3.2.7),  T  =  —  being  the  reduced  electron 

transit  time.  It  is  now  necessary  to  obtain  a  relationship  between  the  field  at  the  entrance 
electrode  Et  and  the  injection  current  l .  As  the  first  step  we  express  Ej  and  l  in  terms 
of  the  reduced  total  transit  time  rd  .  Thus  eliminating  E,  from  (A2.7),  (A2.8)  and  bearing 


in  mind  that  at  the  second  electrode  (collector)  X2  =  1 ,  X‘2  =  1  and  X  =  rd  ,  we  obtain: 


54(-  A 

Td 


(A2.8) 


Similarly  from  (A2.6)  we  obtain 
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(A2.9) 


12 


(rd-l) 


where  we  have  used  (A2.8)  to  eliminate  t .  The  two  equations  (A2.8)  and  (A2.9)  are 
shown  in  Fig.  Al;  they  give  us  the  required  relationship  between  Ej  and  i  in  a 
parametric  form,  the  parameter  being  the  total  transit  time  Td  .  In  the  absence  of  the 
magnetic  field  we  can  go  a  step  further  and  eliminate  ta  altogether  obtaining  a  direct 
relationship  between  Ej  and  l  in  any  of  three  forms: 


27Ef  =(72E1  —  16z)/ 


(A2.10) 


or 


i 


f 

1  + 

V 


(A2.ll) 


or 


^  =  E, 


(A2.12) 


The  above  equation  is  shown  in  Fig.  A2  where  the  electric  field  at  the  entrance  electrode 
E,  is  plotted  as  a  function  of  the  injection  current  l .  As  we  shall  see  later  only  some 
parts  of  the  above  curve  are  of  physical  significance.  In  particular  point  A  indicates  the 
maximum  value  of  the  injection  current  l  =  imax  =  8.0 ,  which  leads  to  a  stable  solution 
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[3].  Similarly  point  B  corresponds  to  an  unstable  solution  for  which  x'B'  =  xB  =  0  or  in 

\ 

=  0  from  (A2.2)  and  (A2.1).  If  we  now  consider  a 

Jb 

curve  which  is  obtained  by  putting  x'  =  0  in  (A2.6),  X"  =  0  in  (A2.7)  and  then  eliminate 
the  transit  time  f ,  we  obtain 


this  case  <j)B 


■  0  and  EB  =  ■ 


d  <f 
dx 


E, 


(A2.13) 


The  curve,  is  shown  dotted  in  Fig.  A2  and  it  crosses  the  original  E,  =  E,  (i)  curve  at 
point  B. 


A2.2  Distance-dependent  solutions 

In  order  to  obtain  solutions  in  terms  of  the  distance  between  the  electrodes  we  substitute 
(A2.2)  directly  in  the  Poisson  equation  by  writing: 


d  V  _  d2u2  _  _  _  Jc  g 

dx2  dx2  P  x-  u 


(A2.14) 


Here  we  have  used  u  =  x'  for  convenience  of  comparison  with  the  results  obtained  for 
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the  more  general  case  of  B  ^  0  .  Multiplying  both  sides  of  (A2.14)  by  — and 

dx 

integrating  we  obtain 


149 


(A2.15) 


1  du2 

2  dx 


=  u^  =  ^  u  +  C, 
dx 


Integrating  once  again  we  obtain 


3^' 


-{(<^u  +  C1)3/2-2C]7^u  +  C1j=^TA/^n  +  ci('y'n-2Ci)  =  x  +  C 

I  3^” 


(A2.16) 


Now  by  definition  u  =  1  when  x  =  0  and  the  C2  constant  can  be  easily  eliminated  by 


writing: 


x,=+^~  y^u+C1(^u-2C1)+T^V^+C,(^-2C1)  (A2-17) 


_2_ 

3dT 


Here  0  <  x_  <  xm  and  xm  <  x+  <  1  ,  where  xm  gives  the  position  of  the  potential 
minimum.  Since  we  are  primarily  interested  in  a  short  circuited  diode,  we  can  go  a  step 
further  and  obtain  a  direct  relationship  between  Cj  and  l  by  noting  that  in  this  case 
U  =  u2  =  1  at  the  exit  electrode  as  well.  We  then  obtain: 


4 1 


'  =  ^=9 


(A2.18) 


As  it  happens  the  curve  l  =  z(Cj)  is  somewhat  unwieldy  and  in  this  case  it  is  much  more 
convenient  to  use  a  new  constants  defined  by  writing 
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Substituting  in  (A2.19)  we  obtain 

i=-rV1"5»(1+25m) 

y/l 


(A2.19) 


(A2.20) 


The  function  /(um)  is  shown  in  Fig.  A3  and  is  a  cubic.  We  find  from  (A2.15)  that 
physically  the  new  constant  Um  represents  the  value  of  velocity  at  the  potential 
minimum,  where  du2/dx  =  &<j>  /dx  =  0  . 

So  far  we  have  presented  well  known  facts  [1,3]  in  a  form  which  is  convenient  for 
comparing  the  results  with  and  without  the  magnetic  field,  the  corresponding  potential 
distributions  being  shown  in  Fig.  A4.  We  find  that  as  i  increases,  the  field  at  the  entrance 
electrode  also  increases  until  it  reaches  a  maximum  value  for  t  =  8 .  This  corresponds  to 
the  lower  branch  of  the  curve  between  the  origin  and  point  A,  Fig.  A2.  If  we  attempt  to 
increase  i  beyond  this  point,  the  potential  minimum  drops  catastrophically  down  to 
zero,  forming  a  virtual  cathode  point  B,  which  instantly  moves  towards  the  entrance 
electrode  (anode)  in  order  to  satisfy  the  Poisson  equation  for  two  space-charge-limited 
diodes  back  to  back;  at  the  same  time  i  splits  into  2i-i7  and  ^  >  only  the  latter  part 
reaching  the  second  electrode  (collector).  By  then  the  continuity  of  t  has  broken  down 
and  (A2.4)  which  is  based  on  it  no  longer  applies. 
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We  can  now  interpret  Fig.  A2  as  follows:  for  0<  i  <4.0  only  a  single  value  of  Ej  satisfies 
the  equations;  for  4.0<  i  <8.0  there  are  two  possible  values  of  E!  but  it  can  be  shown  [3] 
that  the  upper  part  of  the  curve  between  A  and  B  is  unstable,  i.e.  it  can  never  be  reached 
in  practice;  in  fact  the  sudden  formation  of  a  virtual  cathode  and  its  shift  towards  the 
entrance  electrode  gives  rise  to  hysteresis  effects  [1].  The  importance  of  the  dotted-curve 
is  that  all  points  of  the  E,(z)  curve  lying  above  it,  i.e.  all  points  between  B  and  the  origin 
are  of  no  physical  significance,  since  by  then  our  simple  equations  no  longer  represent 
the  system,  a  virtual  cathode  having  formed  between  the  electrodes  and  l  ceasing  to  be 
continuous. 
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